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ABST6ACT 

This thesis is a tutorial on digital signal modulation 
technigues used in satellite communications and includes 
computer simulations of those digital signal modulation 
techniques introduced. The purpose of the thesis is to 
introduce digital signal modulation techniques and through 
the use of computer simulation, generate statistics which 
represent the characteristics of the FFT for the respective 
signal type. Further, an analysis of tne statistics of the 
FFT * s was conducted to determine if there is any relation- 
ship between the components of the FFT of the different 
signals. The statistic used to investigate this possible 
relationship was the F-distribution. The computer simula- 
tion was written and conducted in the FORTEAN programming 
language. A copy of the program, results of the simulations 
and the statistical analysis conducted are included *in the 
appendices. 
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I. INTBODOCTION 



A. BACKGfiOOND 

Since the introduction of the sampling theorem and tae 
matched filter, digital communications techniques have 
developed into a highly proficient discipline. The marriage 
of this discipline with the rapidly expanding space program 
has resulted in communication satellites employing a 
multitude of digital signal modulation techniques- A 
modulation technique is a method of transmitting the 
information contained in a message by varying cr modulating 
the characterist ics of a carrier waveform. These methods 
offer a range of advantages and disadvantages depending on 
the specific characteristics of the modulation technique 
employed. The applications of the various signal modulation 
techniques likewise vary. In order to understand this 
exciting new field it is necessary to look at some of the 
aspects of satellite communications systems in general. 
Then an investigation will be made into the general 
attributes of digital communications and their relationship 
to the satellite system. 

Ihe ability to understand these digital signal 

modulation techniques is the first step in being able to 
intercept, identify and demodulate unknown digital signals. 
These digital signals, transmitted from unknown sources, are 
believed to display frequency characteristics peculiar to 
the modulation technique employed in the encoding process. 
The digital computer offers unique opportunities in 
simulating these modulation techniques, in signal processing 
and in decoding of an intercepted signal- 



B. SPECIFIC GOALS 



It is the specific goal of this thesis to investigate 
the open literature on the topic of digital signal 

modulation techniques in satellite communications- This 
includes a basic understanding of the communications 
satellite system as well as the specific techniques employed 
in signal modulation. Additionally, once a basic 

understanding of these digital signal modulation techniques 
is achieved, computer code in the FORTRAN programming 
language will be developed which simulates these modulated 
signals. The statistics of the time-varying Fast Fourier 
Transforms (FFT) of these simulated signals will be 
investigated. The purpose of this analysis will be to lead 
to fcllow-on research in the area of signal analysis of 
intercepted digital signals whereby they can be classified 
by their FFT as employing a specific digital signal 
modulation technique. 

C. SCOPE OF THE PROJECT 

Althougn an indepth analysis of all the digital signal' 
modulation techniques which will be introduced involves 
considerable higher level mathematics, it is not-- within the 
scope of this project to delve heavily into the mathematics. 
It would be advantageous, however, for the reader to have 
had integral and differential calculus and an introduction 
to statistics. Also a course in electrical engineering may 
prove helpful but is not essential since in actuality it is 
the electrical engineering aspects of satellite 

communications that this tutorial is attempting to present. 
The computer programming will be accomplished utilizing the 
techniques of software engineering and top down modular 
design. Existing blocks of code will be used as modules 
whenever appropriate and available. Again it is the 



ultimate purpose of this project to examine a variety of 

ft 

digital signal modulation techniques and develop computer 
code that simulates common signals used in satellite 
communications that have been produced by one of the many 
digital signal modulation techniques to be investigated- 



II. UHDEHSTAHDING SATELLITE COMMOHICATIOHS 



A. HISTOBI AHO APPLICATIOHS 

Forty years ago, in 1945, Arthur Clarke first envisioned 
the use of space stations placed in geosynchronous orbit for 
communicating to different points on the earth [Ref- 1]. 
Nine years later in 1954, J.R. Pierce of Bell Laboratories 
performed a system analysis on such a communications system 
[Ref- 2]. In 1957, the launch of Sputnik demonstrated the 
feasibility of using a satellite for just such an 
application. However, by 1961 the only satellite 

communications technologies which had been demonstrated were 
the Courier IB satellite, a short-life- active retransmission 
teletype communicati ons satellite in a 1000 km orbit, and 
the Echo I balloon in a 1600 km orbit which demonstrated 
passive reflection of powerful microwave signals from one 
earth station to another. Active microwave communications 
demonstrated in Projects Telstar and Relay were years away 

[Ref. 3]. The first geostationary orbit was achieved by 

Project Syncom in 1963 [Ref. 4 3- 

In 1964, communication organizations from several 
countries joined together to form the international 
organization of INTELSAT (International Telecommunica tions 
Satellite Organization). INTELSAT'S purpose was to develop 
a satellite network which would provide truly global 
communications capabilities. This resulted in the launch in 
1965 of the world's first commercial communications 

satellite, INTELSAT I, also known as "Early Bird". With 
"Early Bird", telecommunications utilizing satellite relay 
were established between the United States and Europe. 
[Ref. 3] 
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The reliability of these satellite systems has improved 
dramatically since INTELSAT I in 1965. Reliability of 
individual links in the system approach 99.99 percent or 
higher. Total system reliability exceeding 99.9 percent is 
common [Ref. 5], making satellite communications more 
reliable than many other modes of communications. In this 
sense, reliability is a measure of the probability that no 
failure will occur in a respective channel or in the system 
during the design life of the satellite. 

In crder for the operability, capability and reliability 
of these systems to have developed at this pace, it was 
necessary for the technologies associated with them to 
develop as rapidly or even more rapidly than the systems 
themselves. Engineering sciences and specifically the 
fields of aerospace and electrical engineering historically 
have required between 7 and 10 years to take a concept from 
operational requirement to full scale operation. This was 
not the case with the concept of communications satellites 
which has seen four generations of INTELSAT satellites 
within a decade's time. [Ref. 3] 

The "Early Bird" satellite weighed approximately 38 
kilograms and possessed limited power and bandwidth capacity 
enabling it to carry only 240 two-way telephone 

conversations. Today's communications satellites represent 
order-of-magnitude improvements in many important operating 
parameters such as power and bandwidth. Additionally, 
increased effective radiated power from the techniques of 
stabilized earth pointing antennas have greatly increased 
the capacity of later generation communications satellites- 
INTELSAT V, the current generation of communications 
satellites, is a three-axis stabilized platform using not 
only the 6/4 GHz frequency band (6 GHz receive/4 GHz 
transmit) as in earlier generation satellites but also a 500 
MHz bandwidth available in the 14/11 GHz band- The 



the 14/11 GHz 



band- 



separation in the receive and transmit frequencies is 
necessary to prevent interference during simultaneous 
operation of the receiver and transmitter. Another factor 
increasing the capacity of present generation satellites is 
the increase in primary power available in the satellite- 
[Bef. 3] 

The technologies which have contributed to the evolution 
of communications satellites come from two primary sources, 
namely technology from the space program of the 1960* s 
supported by NASA and DoD and communications technology due 
largely to commercial and private sector contributions. 
[Ref. 3] 

Electronic devices and components have contributed 
significantly to the rapid development of satellite systems. 
These electronic devices range from something which is now 
considered basic, i. e., the transistor, to devices such as 
Traveling Wave Tube amplifiers, lightweight antennas and 
antenna feed systems. There have been major improvements in 
satellite power sources including more efficient solar cells 
and high storage capacity/lightweight batteries. 
Additionally, two significant developments in electrical 
engineering have made modern day digital signal processing a 
reality. They are the sampling theorem and the matched 
filter. [Ref. 3] 

Communications satellites of the future are likely to 
utilize onboard signal processing. Signal processing 
functions such as signal reshaping, switching and/or 
multiplexing and compression could soon take place on the 
satellite due to the reduction in size and weight of the 
necessary hardware. Operation at over 100 megabits per 
second are envisioned. Onboard signal processing will 
greatly reduce the expense and complexity of earth stations 
thereby making services of a satellite available to a wider 
range of small and geographically disperse users. [Ref. 3] 



Satellite communications hold great promise to provide 
service to a multitude of users over a wide area. 
Applications lie not solely in retransmission but also in 
data collection from that same , large geographic area. 
Military users have definite applications in these areas 
when warning, intelligence and surveillance systems provide 
digital inputs to a central source. Although most of the 
present applications for satellite communications are still 
provided by the government, commercial applications nave 
seen tremendous increases in the last several years. 
[Ref. 3] 

B. ORBITS AND LIHITAIIONS 

Satellite orbits are generally categorized as either 
equatorial (0 degrees inclination), polar (90 degrees 
inclination) or inclined at some angle other than 0 or 90 
degrees relative to the spin axis of the earth as 
illustrated in Figure 2.1 [Ref- 6]. Each satellite orbit 
has a character istic velocity which is dependent on the 
height of the orbit and the orbit* s eccentricity.. 
Eccentricity is a measure of the degree to which the orbit 
approximates a circle, A circle has eccentricity equal to 
zero since the focus is located at the center. See equation 

2.1 and Figure 2,2. 

A communications satellite in elliptical orbit about tne 
earth obeys Kepler's second law of planetary motion. That 
is, a satellite's constant angular momentum about the earth 
means that its areal velocity also remains constant. See 
Figure 2.3. Of particular interest is the satellite 
velocity at apogee (Va) and perigee (Vp) given by equations 

2. 2 and 2.3. 

Note that for a circle, e = 0 and Ra = Rp. Therefore a 
satellite in circular orbit about the earth has an orbital 
velocity as given in equation 2-4. 
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Figure 2.1 Satellite Orbits 

Most communications satellites are placed in a circular 
equatorial orbit where the desire is to stabilize them over 
a fixed point on the surface of the earth called the 
subsatellite point. This type of orbit is called 

geostationary or stationary- Any satellite which has an 
orbital period equal to the period of rotation of the earth 
is called synchronous or geosynchronous. These terras are 
used almost interchangeably in most literature. As 
mentioned for a geosynchronous communications satellite^ the 
orbital period of the satellite^ T, must be equal to the 
period of rotation of the earth- Kepler’s third law of 
planetary motion can be rewritten to yield equation 2-5. 

The period of revolution of the earth for a sidereal day 
is 23 hr 56 min 4 sec. For that given period there is only 
one satellite altitude as expressed by Kepler’s third law. 
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Figure 2.2 Eccentricity 



eccentricity (e) = c/a (egn 2.1) 

* 

4 

c •= distance from focus 
a = semi-major axis 



b = semi-minor axis 
Ra = radius cf- apogee 
Rp = radius cf perigee 



By rearranging terms that altitude is given in equation 2.6. 
Since the orbit is circular, a = Ra = Rp and the height of 
the orbit above the surface of the earth is h = a - Re; or 
35, 804 km. 

One disadvantage of a geosynchronous communications 
satellite is the lack of global coverage. These satellites 
provide excellent coverage of the subsatellite point and 
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Figure 2.3 Kepler's Second Lav of Planetary 



Vd = (k/Ra (1-e) )**. 5 
Vp = (k/Rp (1+e) )** . 5 

* 

4 

k = GMe = 3.98866 x 10**-11 ni3/s2 
G = universal gravitational constant 
G = 6.67 X N m2/kg2 

N = newtons = m kg/s2 

Me = mass of the earth = 5.98 x 10**24 kg 
Vc = (k/Rc)**.5 
Rc = Re + h 

Re = radius of the earth = 6.37 x 10**6 m 
h = satellite orbital altitude 

T = 2 pi a**1. 5 / k**. 5 

T = period of rotation 
a = semi-major axis 



Motion 

(egn 2.2) 
(egn 2.3) 



(eqn 2.4) 
(egn 2-5) 
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a = ,(T / 2 pi) **2/3 (k)**1/3 



(eqn 2-6) 



a = 42,173 km 

laterally to a latitude of about ±80° [Ref. 6 ]• This is 
generally no problem for commercial applications since the 
polar regions do not require a significant degree of access. 
The military, on the other hand, does have an interest in 
communications in the polar region and therefore has several 
communications satellites that have orbits inclined at 
various angles- This type of orbit generally has 

disadvantages of lack of continuous coverage and a much more 
complicated system of ground tracking and receiving 

stations. 

C. FBEQOENCY BAND CONSIDERATIONS 

Although there appears to be an infinite number of 
frequencies available for communications, restrictions do 
exist as to those which are practicable. Limitations on 
available frequency bands for satellite communications are 
due tc the need to select segments of the electromagnetic 
spectrum which reduce noise and interference and are most 
favorable in terms of power efficiency and propagation 
distortions. Trade offs must be made to arrive at the 
optimum frequency for a particular application since single 
frequencies seldom offer the best performance for all 
variables. The problems arise when consideration is given 

to the number of users requiring the same frequency bands 
including terrestrial communications networks. Since the 
problem of interference is global, a worldwide organization 
has been established to assign frequency bands for various 
applications. This organization is called WARC, World 
Administrative Radio Conference [Ref. 63- 



Illustrations of 



the pcrtions of the electromagnetic spectrum in question and 
current eillocations of satellite frequency bands are shown 
in Figures 2.4 and 2-5 [Bef. 6]. 
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Figure 2- 4 Electromagnetic Spectrum 

In general, as long as the hardware and technology will 
support it, higher frequencies are more desirable and more 
in demand because they offer higher tneoretical capacity- 
This is due to the fact that only a percentage of the 
carrier frequency is capable of actually transmitting a 
signal of a given bandwidth. Higher frequencies would also 
experience less interference with existing land and 
satellite systems, A further discussion of the bandwidths 
available as a function of frequency will be included at a 
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Figure 2.5 Satellite Frequencies 



later point under 
techniques. 

As an additional 
to handle more 
electromagnetic spec 
have teen developed 
rreq uency- division 
allocated satellite 
users into specific 
is time-division mul 
frequency band is sh 



the ■ topic of actual digital signal 



note, in order to assist in the ability 
signals in an obviously limited 
trum, several multiple access schemes 
, The most common techniques involve 
multiple access (FDMA) where the 
frequency band is divided among the 
uplink and downlink frequencies. Next 
tiple access (TOMA) where the satellite 
ared by all users by carefully dividing 
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the time any one user has access. Finally there is 
code-division multiple access (CDMA) which involves 
modulating a specific address code which has been 
superimposed on the signal directly onto the carrier. In 
this manner all users share the satellite frequency band and 
only those receiving stations that can demodulate the 
address code can receive the specific signal. [Bef- 6] 

D. WHY DIGITAL MODOLATION 

In general, digital signal modulation techniques are 
superior to analog signal modulation techniques used in 
satellite communications for the following reasons; 

1 . Compatibility with digital computers 

2. Economic advantages 

3. High degree of flexibility 

4. Less susceptible to interference 

5. Quality of signal independent of transmission 
distance and network makeup 

1 . Compa tibilit y with Digita l Com puters 

Clearly one of the most distinctive advantages of 
digital techniques over analog techniques involves computer 
applications. Properly formatted digital signals can be 
used to represent any analog signal (more on this under the 
discussion of the sampling theorem). Once in the digital 
format these signals can be easily manipulated within the 
digital computer. Arithmetic operations can be applied as 
well as logical operations. The signal can he stored 
without alteration or delayed and therefore can be used to 
"simulate" real physical situations. 

The hardware associated with digital circuits is 
free from drift or aging and does not require calibration. 
Additionally digital circuitry is compatible with present 
day integrated circuit technology allowing a standardized 
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building block construction approach. The operating 

• 

characteristics of thes^e systems employing a digital 
computer can be changed by altering the software rather than 
the hardware as is the case fcr analog systems- Finally, 
the use of digital computer technology allows time 
multiplexing. 

2. Flexibil i tv and Econom y 

The flexibility and economy of digital satellite 
communications comes from the fact that more and more 



processing can be done onboard 
downlink to be completely s 
nature makes it possible fo 
reliability through the use 
available to analog systems 
processing or multiplexing is 
signals, simpler and cheap 
stations and terrestrial 
possible. Additionally, ther 
and increased capacity assoc 
[Ref- 7] 

3 . Quality a nd Inte rferen 

The capability of digi 
signal and allow for mul 
processing without degradatio 
digital signal basically 
distance. Multiple hops from 
satellite to satellite are po 
the ncise characteristic of a 
digital systems are capable 
noise ratio of 20 dB to 30 dB 
requiring a much more powerful 



- This allows the uplink and 
eparated- This regenerative 
r low error rates and high 

of digital techniques not 

- Because digital signal 
less costly than for analog 
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III. I MTRODaCTION TO THE FCOfilER TR AHSFOBB . SAMPLING 

THEOREM# 

A DTOCOR RE LATIO N FUNCTION THS ^TCHED FILT ER 

Essential to the understanding of digital signal 
modulation techniques are a few basic tools of the 
electrical engineer and the mathematician. These tools will 
be developed and elaborated on to the extent necessary to 
understand their applicability to the subject of digital 
communications. The description is not meant to be a 
detailed investigation of the respective topics. where 
relevant, the application of the concept being described 
will be mentioned. 

A. FCURIER TRANSFORM 

In mathematical terms, voltages can be expressed as 
functions of time or of frequency. It is more common to see 
voltages represented as a function of time as in equation 

3.1 . 

V (t) = A cos(wt) (eqn 3.1) 

A = amplitude 

w = angular frequency = 2 pi f 
f = frequency = 1/T 
T = period 



It is important to note that both the time and frequency 
functions are representations of voltage and as such may be 
used interchangeably. The Fourier representation [v(t) ] = 
V (f ) is a voltage descriptor in the frequency domain (a 
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function of frequency) while v (t) is a voltage descriptor in 
the time domain (a function of time) . They are different 
but equivalent and either voltage descriptor may be used, 
depending on the concept being explored, to best represent 
the voltage within context. 

The Fourier transform is of principle interest rather 
than the Fourier series since the latter is applicable only 
to periodic voltages. The Fourier transform is applicable 
to voltage pulses, random voltages and other non-periodic 
vol tages, 

DEFINITION : 

CX5 

s^[v(t)] = V (w) = Jv(t)e**-jvt dt (eqn 3.2) 

- oo 

v(t) < > V(wJ 

time domain frequency domain 



Remembering Euler's formula 
e**-jwt = cos(wt) - j sin (wt) 



(eqn 3,3) 



the Fourier transform becomes 



vycv(t) 1 

oo 

V ( w ) = J V 

—oo 



oo 



= V(w) = j'v(t) [cos ( wt) - j 

oo 

(t)cos(wt) ~ jj v(t)sin(wt) 



— oo 



sin (wt) ] 



(eqn 3-U) 
(eqn 3.5) 



Since it will be seen that digital communications deals 
primarily with pulses of finite duration (expressed as 
period, T) , it is worthwhile to examine the Fourier 
transform of a pulse of amplitude A and duration T- 

I A; -T/2 < t < T/2 

1 et V (t) = i 



0 ; el sewhere 
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Figure 3.1 Square Pulse 



Figure 3. 1 is a representation of v (t) or the voltage 
expressed in the time domain. The position of v(t) on the 

A 

t-axis was chosen for convenience of integration but could 
have been situated anywhere on the time line. 



^C»(t) ] = v(w) 

■ f ■ 

V (w) = J e**- ]wt 

-CO 

r 0»e j wt dt 

^2 T/^ 




e**-jv(t dt 
e**-jwt dt 



+ 



(egn 3.6) 
(egn 3.7) 



V(w) = JA e**-jvt dt (e'm 3-8) 

-Vi. 

The integration above may be attacked either head on or 
by substituting cos(wt) - j sin(wt) for e**-jut. The direct 
approach is illustrated due to the relative simplicity of 
the integrand. 

V(w) = Ja e**- jwt dt (egn 3.9) 



T/2 

V(w) = -A/ j w[e**- j w t ] j (egn 3.10) 

-T/2 
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V(w) = - A/ jw[ e**- j wT/2 - e**jwT/2] 



(eqn 3.11) 



V (w) = A/jwC e** jwT/2 - e**-jwT/2] <eqn 3.12) 

By substituting 2pi f = w, the following result is obtained: 

V(f) = A/j2 pi f[e**j2 pi fT/2 - (egn 3.13) 

e**-j2 pi fT/2 ] 

V(f) = A/j2 pi f[e**j pi fl - e**-j pi fT] 

Using Euler’s formula, i-e.^ sin 6 = j9-e**- j9) /2 j ; 

V (f ) = 2jA/j2 pi f[ (e=^*j pi fT - (eqn 3.14) 

e#*- j pi fT) /2 j ] 

V (f ) = A/pi f[ sin( pi f T) ] 

Knowing that sin x / x = sine x 

V (f ) = AT/pi fT[sin (pi f T ) ] (eqn 3.15) 

V (f ) = AT sine (pi fT) 



The sine funetion is eommon in digital eleetronies and 
plots as the produet of sin(pi fT) and 1/(pi fT) as in 
Figure 3.2. In Figure 3.2, V (f ) in the frequeney domain is 
equivalent to v(t) in the time domain. Note that as the 
pulse, T, gets longer, 1/T gets smaller or the first zero 
erossing of the sine funetion oeeurs at a lower and lower 
f re y ueney. 
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1 . Amplit ud e and Phase Spectru m 

Recall that V (w) can be represented as in equation 
3.16. Ihe cosine term is the real part while the sine term 
is the imaginary part. By referencing Figure 3.3, a brief 
review of the complex plane is accomplished and its 
relationship to the amplitude and phase spectrum of a given 
voltage is’ represented. See equations 3.17 through 3.20. 

oo CO 

V(w) = Jv(t)cos(wt) - j Jv (t) sin (wt) (eqn 3.16) 

— CO real ~®^imaginary 



!V(f) I is called the amplitude spectrum of the given 
voltage. The amplitude spectrum can also be calculated by 
using the complex conjugate of the Fourier transform and is 
always positive as shown in equation 3.21. The phase 
spectrum, 0, of the function in question is represented by 
the arctanf imaginary /real ] and is illustrated in Figure 3.5. 
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X^r^ 




Figure 3.3 Phaser Diagram 



r = [ a 2 + b 2 ] 

a + jb = [a2 + b2]**.5 

V (f ) = [ redl2 + imagiaary 2 ]♦*. 5 

I V (f) I = [real2 + imaginary2 ]**. 5 
since |e**j9| = 1 



(eqn 3. 17) 
(egn 3. 18) 
(eqn 3. 1 9) 
(eqn 3. 20) 



1 V (f) I = [ V (f) •V=«'(f )]=«'=«'. 5 (egn 3.21) 

2 - Properties of the Fourier Transform 

Several properties of the Fourier transform are 
useful in the study of digital signals. They are 

represented here without proof and without a great deal of 
det ail. 

a. Linearity 

if v1 (t) <-> V 1 (f) and 

if v2(t) <-> V 2 (f ) then 

■Ji avl (t) + bv2 (t) ] = aV 1 (f) + b/2 (f) 
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Figure 3-5 Phase Spectrum of Square Wave 



b. Time-delay 

^^v(t-to) ] = V (f) e**-jwto 

Note that the amplitude spectrum of the delayed 
version is the same as the amplitude spectrum of the 
undelayed version. It is comforting to note that the 

Fourier transform of a pulse is the same tomorrow as it is 
today, 

c. Scale change 

J^[v(at) ] = VI al V(f/a) 
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d . Frequency translation 

v(t)cos(2 pi fc t) <-> 1/2[V(f+fc) + V(f-£c)] 

V (t) is any voltage 

cos (2 pi fc t) is a carrier wave of frequency fc 

Since understanding of this very important property 
of the Fourier transform is essential to the understanding 
of digital signal modulation it is expanded slightly here. 

As we have seen, the Fourier representation of a 
voltage pulse of amplitude A and duration T is the sine 
function of amplitude AT as illustrated in Figure 3.6. 




Figure 3.6 Square Wave in Time/Sinc Function in Frequency 

The translation is the product of v(t) and in this case 
cos (2 pi fc t) . In the time domain this product only exists 
in the interval between -T/2 and T/2 as in Figure 3.7. The 
significance of this translation and its relationship to the 
bandwidth of the voltage will be discussed further in the 
section dealing with bandwidth- 
e. Differentiation 

d V (t) /dt <-> j w V (f) 

A differentiator could be used as a clock for timing 
but would never be used in the presence of noise since the 
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Figure 3.7 Hultiplication and Translation Diagrams 

result is exaggerated in the presence of high frequencies 
because w = 2 pi f. 
f. Integration 

V (t) dt <-> 1/w V (f) 

An integrator could be used to reduce the effects of 
noise since w = * 2 pi f is in the denominator tending to 
deemphasize the presence of high frequency noise. 

B. THE SAaPLING THEOBEM 

Essential to the understanding of digital c oramunications 
is the sampling theorem which was first introduced by 
Myquist in 1928 [Ref. 9], and later by Shannon in 1948 
[Ref. 10]. The sampling theorem states that any voltage can 
be uniquely represented by appropriately spaced sample 
values of the original voltage. More correctly stated# the 
sampling theorem places limits on the accuracy with which a 
signal can be represented- 

The implication is that an analog signal can be 
represented digitally or by a set of numbers# i-e. # samples. 
A description of the sampling theorem follows. 
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Given any analog voltage, v (t)- as in Figure 3.8, the 
sampldlng theorem says that the entire analog signal is not 
required to accurately represent the voltage but only 
samples of it, call them vs (t) . Figure 3.9 shows samples of 
a representative analog voltage. Samples can be taicen of 
v(t) at every T seconds for a period of seconds. Ihis can 
be accomplished by the use of a voltage clock, call it 
vc(t). The sampling can be viewed graphically as a block 
diagram representing an analog voltage multiplier as shown 
in Figure 3.10. To be of further use in the understanding 
of digital communications we are interested in a frequency 
description of the sample voltage, vs(t). Note that the 
system which describes the obtaining of vs(t) involves a 
voltage multiplication. It was demonstrated in the 

proceeding section that voltage multiplication amounted to 
frequency translation, a property of the Fourier transform. 




Figure 3.8 Representative Analog Voltage 

First it is necessary to find the frequency description 
of the clock, vc (t) , Let vc{t) be a periodic square wave of 
height 1 and duration d as in Figure 3,11. Since vc(t) is 

periodic, it can be shown that the Fourier series 

representing vc(t) is given by equation 3,22, 
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Figure 3.9 Saiples of a Representative Voltage 



r 



•'t 







*• vs(+) 



Figure 3-10 Analog Voltage Multiplier 

Again it is noted that vs (t) , the sample voltage, is the 
product of v(t), the original analog voltage, times vc(t) , 
the clock voltage. In other words vs (t) = v (t) times a 

series cf cosine terms. 

oo 

vc (t) = ao +^an cos (2 pi n fc t) (egn 3.22) 

t\-4 

ao and an are left unevaluated 
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Figure 3.11 k Representative Voltage Clock 



. ^ 

vs(t) = vH) [ao ^an cos (2 pi n fc t) ] (egn 3.23) 

hs I 

Now, assuming any form of the Fourier transform of v (t) , 
as in Figure 3. 12 and since vs(t) is a product and by 
utilizing the frequency translation property of the Fourier 
transform, the following is obtained as a representation of 
the frequency spectrum of vs (t) . See Figure 3.13. 

Remember that the curve highlighted in the box in Figure 
3.13 is the Fourier represen tion of v (t) , the original 
signal. This original signal can now be recovered by 
filtering with an appropriate low pass filter of bandwidth 
equal to or greater than B. This low pass filter would only 
permit the reception of that portion of the signal which 
represents the original voltage. 

The only question left to resolve is how often to take a 
sample. Again referring to the diagram in Figure 3.13, it 
is noted that in order to prevent any overlap of successive 
translations (aliasing) 
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Figure 3.12 Fourier Transform of v(t) 
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Figure 3-13 Fourier Transform of vs (t) 

fo - E » B 
f o » 2B 

or th€ frequency of the clocks fc (the sampling rate = 1/T) 
must be strictly greater than 2 times the largest 
significant Fourier frequency component present- Since 3 
will vary with different v(t)^ the sampling rate necessary 
to uniquely recover that v(t) also varies. 
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In summary on the sampling theorem, it can be said that 
any analog signal can be represented as certain sample 
values spaced the appropriate distance apart. The sampling 
theorem places limits on the accuracy of this 
representation. These sample values can then be transmitted 
from one position to another using analog to digital 
conversion and any one of a variety of modulation 
technigues. It has been shown that the Fourier frequency 
representation of the sample is not equivalent to the 
Fourier frequency representation of the source voltage; 
however, the source voltage can be recovered at the 
receiving end by filtering. A block representation of the 
system is shown in Figure 3.14. 




Figure 3.14 Analog Voltage Multiplier and Low Pass Filter 



C. TEE AOTOCOERELATICH FONCTION 

The above concept certainly seems simple enough. If a 
signal is present and it is desired to transmit it from 
point A to point B, all that has to be done is to take 
appropriately spaced samples of the signal, encode them. 
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somehow transmit them to point B, filter the received signal 
and the uniqueness of the original voltage has been 
reproduced. Unfortunately it isn’t quite that easy. The 
reason it is not, is due to the presence of noise. The 
sources of noise will not be discussed here, however, the 
effects of noise in general terms and how a faint signal can 
be recovered in the presence of noise will be discussed. 

The time varying descriptor of voltage and the frequency 
varying descriptor of voltage have been introduced. Both 
are precise mathematical descriptions of something (voltage) 
that is deterministic. That is, it can be described in 
mathenatical or graphical terms during any period of time. 
Such is not the case for noise corrupted voltages which are 
entirely random. Therefore other descriptors of voltages 
must be employed in the presence of noise. They are partial 
descriptors of voltage since a random signal cannot be 
described with precision. These partial descriptors are the 

1. Autocorrelation function and the 

2. Probability density function (p.d.f.) 

It is important here to note the difference between the 
source voltage at the receiver and the sample voltage. The 
source voltage at the receiver, vsr(t), is the digitally 
converted sample voltage, Vs(t), modulated onto a carrier by 
one of the digital modulation techniques yet to be 
discussed. understanding the autocorrelation function is 
the basis for understanding how a decision is made that 
vsr (t ) is present at the receiver in the presence of noise. 
It should be remembered that the exact form of the carrier 
is known at both the transmitter and receiver. 

First of all, an additive noise model is assumed where 
the signal at the receiver, vr (t) , is equal to the signal 
which is being transmitted, vsr(t), (i-e., the digital 

samples of v(t) modulated onto the carrier), plus random 

noise, vn(t) 



40 



vr (t) = vsr (t) + vn(t) 



(egn 3.24) 



vr (t) = voltage at the receiver 

vsr (t) = signal voltage at the receiver 

vn (t) = noise voltage at the receiver 



If it is realized that most receivers operate on the 
principal of detection of DC voltage, the problem becomes 
one of rectification of the received signal and 
determination if the DC voltage, characteristic of the 
transmitted signal is present. 

A common type of rectification of an analog voltage 
involves squaring the input waveform. If the given signal 
at the receiver is vr (t) as in equation 3.24 above, squaring 
the waveform introduces the square of not only the desired 
signal but also the square of the noise term. If instead of 
squaring vr(t) and introducing or at least not eliminating 
the noise, vr (t) is multiplied by vsr(t), the results in 
equation 3.25 are obtained. 



vr(t)*vsr(t) = vsr**2(t) + vsr(t)vn(t) 



(egn 3.25) 



By averaging, all that remains is vsr **2 ( t) since the 



average of vsr(t)vn(t) is zero because vh ( t') is random and 
the average of any random voltage is zero. The DC component 
of vr (t) «vsr (t) is vsr**2 ( t) - Any remaining AC component 
can be removed by a low pass filter. Graphically in block 
diagram form this receiver looks like Figure 3-15 
illustrated below. If vsr**2 (t) > 0 then vsr(t) is present 
in the signal vr(t). If vr (t) is pure noise or some other 
signal is present exclusively, then the output of the 
receiver will be 0, 
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Figure 3.15 Rectification through an AVH and IFF 

In practice, however, the exact waveform of vsr (t) may 
not be known or it may be a time delayed or distorted 
version of the original signal when it arrives at the 
receiver. If the distortion or delay is significant enough 
there will be little agreement or ’’correlation" in the 
receiver. 

The solution is to multiply the received signal by a 
series of time delayed approximations of the transmitted 
signal. The signal is present whenever the output of the 
series of voltage multipliers or correlators exceeds a 
certain threshold approaching vsr**2(t). The concept is 
illustrated in Figure 3,16, 

For the concept illustrated in Figure 3.16 to work, the 
average value of vr ( t ) «vsr ( t-p) must have a DC component. 

vr (t) «vsr (t-p) = vsr (t) vsr (t-p) + (eqn 3,26) 

vn (t) vsr (t-p) 



vr (t) «vsr (t-p) = vsr (t) vsr (t-p) + 



4 2 



vn (t) vsr (t-p) 



(egn 3. 27) 




Figure 3.16 Correlation Device 

But the average of vn (t) vsr (t-p) = 0 since vn (t) is random and 
the results of equation 3.28 are obtained. 



vr (t) *vsr (t-p) = vsr (t) vsr (t-p) (egn 3.28) 

What is of interest is the average value of 
vsr (t) vsr (t-p) . The autocorrelation function (ACF) will be 
defined as that average value of vsr (t) vsr (t-p) . The 
mathematical representation of the autocorrelation function 
for a continuous voltage is represented in equation 3.29 
while the autocorrelation function for a pulse voltage is 
presented as equation 3.30. 

r 

ACF = Rvv(p) = 1/2T j v(t)v(t-p)dt (eqn 3,29) 

-T 

oo 

ACF = Cvv(p) = J v( t) V (t-p) dt (eqn 3-30) 

- CO 

The autocorrelation function is a measure of the degree 
to which two identical signals which are corrupted in some 
manner are alike. The ACF of two signals which are 
identical, not distorted in any way and occurring at exactly 
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On the 



the same time, i.e. , p = 0 has a maximum value, 
other hand there will be very little correlation between two 
identical time variant signals when the time difference 
between them is great. In this case the autocorrelation 
function is near zero. 

A similar concept is used within the context of signal 
comparison. This concept is the crosscorrelation which is a 
measure of the degree to which 2 different signals are 
alike. When the crosscorrelation between voltage vsr and 
random noise voltage vn is 0 there is no relationship or 
correlation. 

Again let us assume that the signal at the receiver is a 
noise disrupted version of the signal originally 
transmitted. 

vr(t) = vsr (t) + vn(t) (eqn 3.31) 

In the receiver , the time delayed version of the signal is 
applied to the incoming signal and the average is formed as 
before and shown in Figure 3.17. 
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Figure 3.17 Voltage Correlator 
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C vsr(t) + vn (t) ]•[ vsr (t-p) + vn (t-p) ] = 



(eqn 3. 32) 



vsr (t) vsr (t-p) + vn (t) vsr (t-p) + 

vsr (t) vn (t-p) + vn (t) vn (t-p) = 



Rvsr vsr(p) + Rvn vsr (p) + Rvsr vn(p) + Rvn vn (p) 



Rvn vsr(p) = 0 



Rvsr vn(p) = 0; because average vn = 0 



Therefore the result is simply Rvsr = vsr(p) + Rvn vn(p) 
or the autocorrelation function of the desired signal plus 
the autocorrelation function of the noise. Assuming the 
shape of the autocorrelation function of both components 
(both the original signal and the noise which vary with p, 
time) is known, what is done in practice is to vary the value 
of p until the ratio cf Rvsr vsr(p) / Rvn vn (p) is a maximum 
or the autocorrelation function of the signal is a maximum 
while the autocorrelation of the noise is minimum and tne 
signal is recovered. 

D. THE MATCHED FILTER 

The matched filter is a linear filter which has the 
characteristic response desired for optimum reception of the 
desired signal. In other words, we desire the response of 
the system to the linear filter to be in some way 
proportional to the autocorrelation function of the desired 
signal and in no way related to the autocorrelation function 
of the noise. This system could be illustrated as in Figure 
3.18. 
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Figure 3.18 Matched Filter Block Representation 

The question is what should the makeup or h (t) ^ (i.e., 

the response of the linear filter) be? The output of a 
linear system is a convolution so the output of the linear 
filter, h (t) to an input vsr (t) would be vsr(t) * h (t) , 
where * denotes convolution - 

In the previous section it was stated that the 
autocorrelation function, Rvsr vsr(p), a function of p, is 
given by 

CO 

Rvsr vsr (p) = Jvsr (t) vsr (t-p) dt (eqn 3,33) 

- oo 

By simply changing variables, equation 3.33 can be rewritten 
as a function of time, t. 

oo 

Rvsr vsr (t) = vsr (p) vsr (p-t) dp (egn 3.34) 

- oo 

From convolution theory it is known that the response of a 
linear filter as above, vsr(t) * h (t) can be expressed 
mathematically as in equation 3.35, 

« 

oo 

vsr (t) * h(t) = Jvsr (p) h (t-p) dp (eqn 3,35) 

— oo 
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As stated earlier, the desired output of the linear 
filter is the autocorrelation function Evsr vsr(t). Note 
the similarities between the last two equations. If in 
equation 3.35 a simple substitution of h (t) = vsr(-t) is 

performed, the desired results are obtained. 



= R vsr vsr (t) 

Therefore, for every signal that is desired to be 
recovered, the matched filter will do the job very nicely if 
the response of that filter, h (t) is equal to the inverse of 
the signal that is desired to be detected. 



oo 



vsr (t) ♦ h (t) 




dp 



(eqn 3.36) 



-OO 



I?. MliOG TO DIGITAL CONVEBSION 



Key to the understanding of digital communications is 
how the analog information is converted to digital 

information. As illustrated earlier, each analog voltage or 
signal can be represented by appropriately spaced sample 
values. These sample values are still analog in that they 
can take on any value in the range of the original analog 
signal. What is desired is to change this infinite set of 

decimal numbers to a finite set of decimal numbers. This is 
called analog to digital conversion. Common types of A to D 
conversion include Pulse Code Modulation (PCM) , Differential 
Pulse Code Modulation (DPCM) , Delta Modulation (DM) , Pulse 
Amplitude Modulation (PAM), Pulse Duration Modulation 

(PDM),and Pulse Position Modulation (PPM). PCM, DPCM, and 
DM will be examined in this chapter because of their 
widespread usage and easy application to digital technology- 

A. POLSE CODE MODOLATION (PCM) 

PCM is the most widely used A to D conversion technique. 
It involves assigning the sample value obtained in sampling 
to one- of several quantization levels within the range of 
the voltage sampled and then representing those quantization 
levels by various binary code words. 

For example, the highest significant frequency component 
in human speech is 3300 Hz. In order to capture accurately 
the signal produced in speech, one would require a sampling 
rate greater than or equal to 2B or 2(3300) = 6600 

samples/sec. The telephone company uses 8000 samples/sec. 

B = 3300 Hz 

f > 2B = 6600 samples/sec 

use f = 8000 samples/sec to ensure no aliasing 
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If a signal is sampled at the rate of 8000 samples/sec, then 
in the transmission of those samples, assuming the samples 
are transmitted immediately and not delayed, the 

transmission rate must be 
Ts = 1/f 

Ts = 1/8000 samples/sec 
Ts = 125 microsec/sample 

The number of quantization levels employed to represent 
the various analog sample values is arbitrary. The greater 
the number of levels, the more accurate the reconstruction 
of the original signal but the more rapid the data 

transmission rate must be (in all cases there will be some 
error present after recovery). Consider again the example 
of voice as illustrated in Figure 4.1. 




Assume it is desired to represent each sample value (decimal 
number) with an 8 bit binary code word. Then the number of 
quantization levels is given by equation 4.1. 

The spacing between quantization levels is the range of 
voltages to be represented divided by the number of 
quantization levels. The analog to digital conversion takes 
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Figure 4. 2 Analog Nature of Samples 
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bit rate = (n bits/ sample) (1 sample/T sec) (eqn 4.2) 

bit rate = (8 bits/sample) (1 sample/125 microsec) 

= 6 4 k bp s 

bit duration = 1/bit rate (eqn 4.3) 

= 1/64 kbps = 15.6 microsec/bit 



( i^signals/channel) (bit rate/signal) = (egn 4.4) 

bit rate/channel 

In other words , if 24 voice signals are to be 
multiplexed onto a single channel, the bit rate must be 

bit rate/channel > (eqn 4.5) 

(24 signals/Channel) (64 kbps /signal) 

= 153 6 Mbps 

As can be seen, the data transmission rate increases 
significantly, necessitating better and better hardware and 
a wider and wider bandwidth. Advantages should be readily 
apparent for A to D conversion techniques or modulation 
techniques that reduce the bit rate required. 

Recovery of the analog signal is by a process called 
Digital to Analog Conversion. All that will be said about D 
to A conversion is that it is the inverse process of A to D 
conversion and that a slight error is always introduced 
during the process due to the discrete nature of the 
quantization process during A to D and D to A conversion. 

B. DIFFERENTIAL POLSE CODE MODULATION (DPCM) 

In BPCM, what is converted to binary numbers is not the 
quantization level (decimal number) which each sample value 



is represented by but the successive differences between 
quantization levels. The idea is that the range of maximum 
difference values will be smaller than the range of actual 
sample values. It is therefore possible to represent that 
range of delta values with fewer bits/binary word and 
therefore fewer quantization levels and ultimately a lower 
bit rate. [Bef. 11] 

In actual practice, in a DPCM conversion technique, it 
is a statistical estimate of each successive sample value 
which is subtracted from the actual value that is converted 
to binary code words. The result is the same in that the 
range of amplitudes is reduced and therefore fewer bits/word 
are required to represent the sample thereby lowering the 
data transmission rate required to transmit the signal. 
[Ref. 11 ] 

C. DELTA lODOLATION (DH) 

Delta modulation is a form of DPCM where successive 
quantization levels of the output differ by only 1 bit. 
That is to say that any successive quantization level can be 
represented by varying only one bit of successive output 
binary code words. In other words a type of gray code is 
employed. This A to D technique is implemented through the 
use of a DM coder or linear delta modulator. This DM coder 
approximates a given input signal with a series of linear 
segments of uniform slope. A comparison is made between the 
value of this approximation and the input signal at each 
sample increment. The sign of this difference value is what 
is encoded and is used to increment the DM coder in the 
direction of the input signal. By using the differential 
sign value of the input signal and the incremented 
approximation from the DM coder the linear approximation 
from the linear delta modulator is said to "track" the input 
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signal. (fief. 8] Slope overload of this type of modulation 
technique occurs when the slope of the incoming signal 
exceeds the ability of the DM system to follow the source at 
the sampling rate being utilized. (Ref. 11 ]• 
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V. DIGITAL SIGNAii MODOL AT ION TECHHIQD ES 



Digital signal modulation techniques are the methods of 
encoding information for transmission utilizing digital 
technology. Factors affecting digital modulation include, 
but are not limited to, the physics of the method, hardware 
requirements, bandwidth considerations, power requirements, 
data transmission rates and error probabilities. It is 
these aspects which will be explored in the following 
sections. 

A. DIGITAL HODOLATION PORHATS 

Modulation is the technique by which the characteristics 
of one waveform (called the carrier) are varied or modified 
by the characteristics of another (called the source) . The 
carrier waveform of interest is the sinusoid. It should be 

4 

obvious that the attributes of a particular sinusoid . that 
differentiate it from every other sinusoid are its 

amplitude, phase and frequency. It follows that the 

characteristics of the carrier waveform to be varied by the 
source are its amplitude, phase, frequency or any 

combination of the three. This gives rise to the broad 
general formats of Amplitude Shift Keying (ASK) , Phase Shift 
Keying (PSK) and Frequency Shift Keying (FSK) . All other 
digital modulation techniques are variations of or 

combinations of these basic formats. Other factors involved 
in the description of the digital modulation technique being 
employed include the number of bits being encoded at one 
time, the employment of error correction techniques and the 
baseband (source) waveform. 
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when each bit of the baseband waveform is individually 
encoded by any of the techniques previously mentioned (ASK, 
PSK, FSK ) , the technique being utilized is referred to as 
binary encoding. When more than 1 bit of the source code is 
modulated onto the carrier at one time it is called block 
encoding. Block encoding allows for one of m = 2**k 

waveforms where k = number of bits. 

This paper will not go deeply’ into signal detection or 



demodulation techniques, howe 
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satellite communications, most 
be considered "off the sh 
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communications include the s 
multiplexer and transmitter, 
above hardware requirements n 
digital formatting, however 
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ver, a basic understanding of 
and again further defines the 
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d to as either coherent or 
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Figure 5. 1 Basic Hardware Component Block Diagram 

1 . Sampler 

The sampler is the device or component which samples 
or extracts those single character istic values of a 
naturally occurring analog signal which are ultimately 
encoded into a digital word by the encoder. The frequency 
at which this sampler must operate was derived in the 
discussion of the Sampling Theorem. It was determined that 
the sampler must operate at a frequency greater than 2 times 
the highest significant Fourier frequency component of the 
source voltage. Modern solid state devices capable of 

taking thousands of samples/sec are available at modest 
costs. 

2 . Encoder 

Often referred to as the A to D converter (analog to 
digital converter), the encoder transforms the samples of 
the analog signal derived from the sampler into a digital 
format through one of the analog to digital conversion 

techniques described in the chapter on A to D conversion. 

These digital bits can be stored for later use, coded. 
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delayed or used immediately either individually or in groups 
to modify one of the characteristic qualities of the carrier 
waveform. When sample values of the analog source signal 
are converted into digital bits of information, they are 
simply that, digital bits of information. The carrier can 
only be modified to represent the source information by 
interaction with another voltage. Therefore the value of 
the digital bit (binary) , either 0 or 1, is used to generate 
a baseband waveform or voltage which does the actual 
modulation of the carrier waveform. 

a. Common Baseband Waveforms 

It should be obvious that since it is desired to 
represent binary code words with a representative baseband 
waveform what is required is two levels of voltage. There 
are two basic logic schemes for representation of the 
baseband waveform. They are bipolar or unipolar logic. 
Bipolar logic involves representing the 0*s or 1*s of the 
binary codeword as either +V or -V as illustrated in Figure 
5.2. 
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The 0 or the 1 of the binary code word can be represented 
either as -V or -V, +V respectively. This is a matter 
of convention but must be clearly understood in the various 
component designs in order to ensure compatibility between 
parts of the system. 

Unipolar logic utilizes a voltage to represent 
either of the possible binary digits and the absence of 
voltage to represent the other- Two common conventions of 
unipolar logic are represented in Figures 5.3 and 5.4. 
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Figure 5-3 Unipolar Positive Logic 

Again, which form of unipolar logic might be employed in the 
A to D converter is a matter of convention. 

Variations of these two basic encoding formats 
have the advantages of ease of generation and improved 
decoding and clock recoverability. Two variations commonly 
used in satellite communications are the Non Return to Zero 
(NEZ) and the Manchester waveforms. The NRZ waveforms is 
simply the voltage representation which corresponds to the 
stream of bits represented in the logic scheme chosen. 
There is no transition as long as the same bit is present. 
Choosing bipolar logic as an example, a NRZ waveform can be 
described schematically as in Figure 5.5. 
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Figure 5-5 Non Bet urn to Zero Waveform 

The Manchester waveform of the same source 
voltage is represented with a transition at the midpoint of 
the bit duration from either +V to -V or -V to +V. This 
form cf coding has the advantage of clock synchronization in 
all cases of digital encoding. Schematically, the 
Manchester code can be represented as in Figure 5.6. 
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Figure 5.6 Kan Chester Vaveform 
3 . Modulat or 

The modulator in a digital communications system is 
more commonly referred to as a modem 

(modulator/demodulator) . This device does the actual 
transformation of the digital information into a waveform 
that can be transmitted from one point to another. , The 
modulation as mentioned earlier can modify either the 
amplitude, phase or frequency of the carrier wave. In 
simple terms, the digital information derived from the 
analog source is mapped into the carrier wave. It is this 
mapping that determines the power and bandwidth 
characteristics of the modem. 

^ • Mult i ple xer 

Multiplexing is a common method for increasing the 
utility of a given communications channel. Variations in 
multiplexing techniques give rise to the multiple access 
techniques employed in satellite communications. 
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a. Frequency-Division Multiplexing (FDM) 

Frequency-Division Multiplexing is a technique 
wherety the capacity of a communications channel is 
increased by adding one signal to another. These signals 
occupy discrete nonoverlapping frequency bands. A specific 
signal is recovered by use of a band pass filter for the 
frequency band desired. [Ref. 7] 

b. Time-Division Multiplexing (TDM) 

Time-Division Multiplexing is a technique used 
to increase the capacity of a given communication channel by 
adding signals in time. That is, specific signals or 
portions of signals are allocated specific time slots or 
portions of the carrier wave. These time allocations cannot 
overlap and the signal is recovered through proper 
synchronization of the recovery hardware with the 
multiplexer. [Ref- 7] 

c. Code- Division Multiplexing (CDM) 

Code- Division Multiplexing is a technique for 
increasing the capacity of a given communica tions channel 
through the assignment of a charac teristic code to the 
digital signal before it is multiplexed in the time or 
frequency domain. Since the code used in multiplexing is 
known at the demultiplexer, recovery of the digital source 
code could be accomplished through the use of a matched 
filter or correlator. [Ref. 7] 

C. BANDWIDTH 

For purposes of discussion and for uniformity, when 
examining the various digital signal modulation techniques 
of interest, the bandwidth will be defined as the highest 
significant Fourier frequency component of the signal in 



question. For digital communications, it can be shown that 
by the above definition, the bandwidth of the baseband 
waveform is 1/Tb, where Tb is the bit duration. For 
example : 

1. Voice sampled at the rate of 8000 samples/sec 

2. Each sample represented by an 8 bit code word 

3. Rb = bit rate = (8000 samples/sec) (8 bits/sample) 

Eb = 64 kbps 

4. lb = bit duration = 1/Rb 
Tb = 15.6 microsec 

The Fcurier transform of the baseband waveform of duration 
15.6 microsec is represented in Figure 5.7. 
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Figure 5-7 Fourier Transform of Baseband Waveform 



It should be readily apparent that as the bit rate of the 
digital baseband waveform increases, so does the bandwidth. 
It is desireable to keep the bandwidth of the signal to a 
minimum due to limitations on availability of the 

electromagnetic spectrum. In addition, interference with 
signals which occupy adjacent frequency bands is reduced and 
propogation limitations of the medium and hardware 

limitations such as self- ind uctance and capacitance are 
minimized. 



D. SPECIFIC TECHHIQDES 



Now with some of the basics of digital signals already 
introduced it is possible to investigate some of the more 
significant modulation formats employed in digital satellite 
communications. The organization of this section will be 

around the three general types of modulation, i-e. , phase, 
amplitude and frequency. Within each broad category, 
several techniques will be described under the subareas of 
binary encoding, where the information is modulated onto the 
carrier tit by bit, and block encoding, where groups of bits 
do the modulation. Variations of these areas will be 
pointed out. The charact eristics of the modulation 

technique will be introduced including, where practical, an 
analytical description, phasor representation, error 
probability, spectral efficiency, advantages and 
disadvantages. 

1. Ph ase Sh ift Ke^inq (PSK) 

Phase shift keying is a generic form of signal 
modulation where the baseband waveform is used to modify or 
change the phase of the carrier. Phase, as stated earlier 
is one of the characteristics which distinguish one sinusoid 
from another. In general, phase shift keying offers the 
advantages of power and bandwidth efficiency. 

a. Binary Phase Shift Keying (3PSK) 

As the name suggests. Binary Phase Shift Keying 
(BPSK) results in a modulated carrier waveform consisting of 
two phases of the same sinusoid. The baseband waveform is 
used to change the phase of the carrier bit by bit from one 
phase to the other. In the case of bipolar logic, the bit 
stream consists of ±1. The general waveform of the carrier 
can be represented as a cosine wave of amplitude A, angular 
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frequency w and initial phase angle delta as in equation 

5.1 . 



vc(t) = A cos(wt + delta) (egn 5.1) 

Since two phases are represented in vc (t) , it is 
customary to let them differ by pi radians- This can be 
represented by the modification of equation 5.1 to account 
for a pi phase shift depending on whether the bit, v (t) to 
be modulated is either a +1 or -1 as shown in equation 5.2. 

VC (t) = A cos(wt + delta + pi/2 ( 1-v (t) ) ) (eqn 5.2) 

Equation 5.2 can be expanded to result in a simplified form 
showing that the modulated carrier waveform is simply a 
product of the baseband waveform and the carrier. See 
equation 5.3 where v (t) represents the baseband waveform- 

vc (t) = A V (t) cos (wt + delta) (eqn 5.3) 

A phasor diagram can be constructed to represent 
the modulated carrier waveform in BPSK and is shown in 
Figure 5-8. As can be seen by the phasor diagram in Figure 
5.8, the phase of the modulated waveform depends on the 
value of the baseband waveform and differs by pi depending 
on whether the value of the modulated bit is ±1. 

Remembering the section on the properties of the 
Fourier transform, it is possible to analyze the spectral 
and power efficiency of BPSK. Recall BPSK can be created by 
multiplication. Multiplication of two voltages is covered 
by the frequency translation property of the Fourier 
transform. Frequency translation results in double the 
bandwidth or spectral requirement of the translated waveform 
and half the power. In the case of a NRZ baseband waveform. 
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Figure 5.8 Phaser Diagram of BPSK 

the bandwidth of the non- translated pulse, determined from 
the largest significant Fourier frequency component is 1/Tb, 
where Tb is the bit duration. The bandwidth of the 
translated voltage (BPSK carrier) is 2/Tb. The power in a 
sinusoid is A^/2 and none of this power capacity is lost in 
BPSK since only the phase of the original sinusoid is 
changed. For a Manchester baseband, waveform, the 
non- translated signal has bandwidth 2/Tb and the translated 
wave has bandwidth 4/Tb. 

The advantage of BPSK is that it is relatively 
efficient in bandwidth utilization and power capacity. It 
does have the disadvantage of fairly large sidelobe 
components which contribute to interference with adjacent 
frequencies. This usually necessitates some type of 
filtering before the signal can be transmitted. 

b. Differential Binary Phase Shift Keying (DBPSK) 

The recovery of a Binary Phase Shift Keyed 
signal requires the use of a phase coherent reference signal 
as described in the previous sections on autocorrelation and 
matched filters. This may not always be practical or 
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possitle so a technique called Differential Binary Phase 
Shift Keying has be developed. DBPSK uses the same carrier 
type as BPSK and offers many of the same advantages and 
disadvantages. As stated, the primary difference is the 
lack of necessity for a phase coherent reference. 

The baseband waveform, v(t), for DBPSK is 
generated by comparing the phase of the next bit to be 
transmitted to that of the previously transmitted bit. For 
example, if a +1 is encoded onto the carrier as v (t) , it 
represents a particular phase cf the carrier. If the next 
bit to be sent is also a +1, there is no phase change, i.e. , 
the bit remains the same over two successive bit durations. 
If, however, the second bit to be sent is a -1, this 
represents a pi phase change and is decoded as a bit change 
over successive bit durations. This can be expressed 

mathematically as in equation 5.4 where bk represents the 
bit which will be encoded on the carrier and bk- 1 represents 
the previous bit already encoded and ak represents the 
present bit. 

bk = bk-1 • ak (egn 5.4) 

As can. be seen when bk = -1, there is a phase 
change and the bit changes from one bit to the next 

depending on the logic type being employed. If bk = +1, 
this represents the phase remaining constant over two 

successive bit durations, i.e., no bit change. See Figure 
5.9 [Bef. 8]. Recovery is effected by correlation of the 
received bit to the previously received bit. 

Differential Binary Phase Shift Keying (DBPSK) 
differs from Differentially Encoded Binary Phase Shift 
Keying (DEBPSK) in the recovery step only. DBPSK takes 
advantage of the modulation of phase shift in the recovery 
process while DEBPSK still utilizes a phase coherent carrier 



66 



1 0 1 1 0 1 0 0 0 1 1 

Message 
sequence ’bj^| 



Differentiallv encoded 
sequence ‘ 6 ^ 

Reference digit 

Level-shitted 

sequence 

-1 — — I 1—1 — I 

T ransmitted OrrrrrrOOrrOTTTrrrO 

phase 





Figure 5.9 Differential Binary Phase Shift Keying 

A 

signal in demodulation. Both cases use the same scheme for 
differential encoding. 

c. Orthogonal Binary Phase Shift Keying (Orthogonal 
BPSK) 

In Binary Phase Shift Keying there is always a 
probability, with the interjection of noise on the 

transmitted signal that an individual bit will be decoded 
incorrectly from that which was actually transmitted. In 
some environments, this probability is so high that special 
encoding techniques must be employed in order to reduce the 
probability that a bit or a binary code word will be decoded 
in error. One such technique is Orthogonal BPSK. 

Orthogonal BPSK is a type of group encoding scheme where 
groups or blocks of bits from the original signal are 
assigned redundant bits according to a predetermined 
assignment routine. This new expanded sequence of bits is 
then transmitted as before in BPSK. 
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For example, assume that it has been determined 
that the respective values of the original signal will be 
represented by a k bit binary code word. It has been 
demonstrated that there exists 2**k binary code words, k 
bits long in this type of arrangement. With Orthogonal BPSK 
there also exists 2**k orthogonal sequences of bits (called 
channel symbols, bauds or chips) ffief. 6]^ which represent 
the k binary code words. However, the number of bits or 
chips in the sequence is increased to n = 2**k, thereby 

providing the desired redundancy. Each of the 

sequences of n chips is constructed to be orthogonal to 
every other sequence in order to assure maximum separation. 
See Figure 5.10 [Ref. 6]. 
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Figure 5.10 Orthogonal Binary Phase Shift Keying 

Decoding is also done by blocks through a bank 
or 2**k correlators. Since the decoding is done in blocks, 
the probability that the entire transmitted code word will 
be incorrectly decoded is reduced even though individual 
bits are incorrectly decoded. This is the major advantage 
of Orthogonal BPSK. Figure 5.11 [Ref. 6], shows the 

reduction in error probability with Orthogonal BPSK 
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Figure 5.11 Probability of Error for Orthogonal BPSK 

The major disadvantage to Orthogonal BPSK is 
that, it results in a reduced data bit rate at a given 
transmission rate or the necessity for a higher transmission 
rate if the same bit rate is to be maintained. For example, 
suppose the bit rate of a PSK signal is R bits/sec. If each 
of the k bits in the PSK signal are represented by 2**k 
chips as is the case in Orthogonal BPSK modulation, then the 



bit rate must be (2**k/k)R if the same data rate is to be 
maintained. Additionally, since the necessary bit rate 
increases, the bit duration decreases. It was shown in the 
section on bandwidth that as bit duration decreases as the 
result of a higher bit rate, an increase in bandwidth is the 
result. Therefore, the tradeoff in decreased transmission 
error comes at the expense of bandwidth and bit rate. 
[Ref. 6] 



d. Quadrature Phase Shift Keying (QPSK) 

Quadrature Phase Shift Keying derives its name 
from a four-phase modulation of the carrier waveform. These 
modulations are achieved by simultaneously modulating the 
inputs from two bit streams onto a single carrier. The two 
bit sequences could be from two separate sources or 
successive bits from a single source. To take advantage of 
the orthogonality of the sine and cosine functions, QPSK can 
be viewed as in equation 5.5 as the sum of a BPSK modulated 
sine and cosine. 

VC (t) = A v1 (t) cos(wt + delta) + (egn 5.5) 

A v2(t) sin (wt + delta) 



The signal can also be represented as in 
equation 5.6 where the phase angle, theta, of the modulated 
carrier is shown in equation 5.7. 

vc(t) = 2**. 5 A cos(wt + theta delta) (egn 5.6) 

theta = arctan (v1 (t)/v2 ( t) ) (eqn 5.7) 

The phases possible as a result of substitution of the 
possible bits in bipolar logic in equation 5.7 are ±45 
degrees and ±135 degrees. See Figure 5.12. 
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Figare 5.12 Hodulation and Phases of QPSK 

Since two bits are being encoded at one time for 
transmission, the bit rate in QPSK is twice that of BPSK. 

A 

The BPSK bit rate is equal to the QPSK symbol rate, (symbol 

rate being the time the quadrature bits exist on the 
carrier) . It is the symbol rate in QPSK which determines 
the bandwidth and since symbol rate is equal to bit rate in 
BPSK, the bandwidth of the two modulation techniques is 
identical. The great advantage of QPSK lies in the 
bandwidth utilization. Two tits of information are 

transmitted thereby effectively doubling the bit rate at no 
additional cost in bandwidth. A phasor representation of 
QPSK is illustrated in Figure 5- 13. Each quadrature 

component of the modulated signal contains half the power of 
the tctal carrier or (A^/2) /2 = A2/4. 

e. Offset Quadrature Phase Shift Keying (OQPSK) 

In Quadrature Phase Shift Keying, both bits of 
the two baseband waveforms, which are to be modulated onto 
the carrier, change at the beginning of the respective 
symbol time. This allows the phase of the carrier to change 
up to 180 degrees for each symbol as illustrated in the 
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Figure 5.13 Phasor Diagram of QPSK 

phasor diagram in Figure 5.13. The phase change of the 
modulated signal can be limited to 90 degrees through a 
modulation technique called Offset Quadrature Phase Shift 
Keying. This is accomplished by delaying the application of 
the second baseband channel for Ts/2 seconds, vhere Ts is 
symbol duration as illustrated in Figure 5-14. Modulation 
is then accomplished as before in QPSK. 




Figure 5-14 Offset Quadrature Phase Shift Keying 
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The result is a QESK modulated signal in which 
the maximum phase shift between successive bits is 90 
degrees. The bandwidth and power spectra of an OQPSK 
modulated signal are the same as that of a QPSK modulated 
signal. The advantage of OQPSK is in the limit which is 
placed on phase shift during encoding which simplifies the 
encoding hardware. Additionally OQPSK has spectral and 
interference advantages during decoding [Bef. 6 3- 
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Figure 5,15 Phasor Diagram of OQPSK 



f. Multiple Phase Shift Keying (MPSK) 



Another form of digital modulation is known as 
Multiple or M-ary Phase Shift Keying (MPSK) . Once again, as 
in Orthogonal Binary Phase Shift Keying, k data bits are 
transmitted. The k bits represent the word length and there 
are m = 2**k different binary words, k bits long, for this 
type of modulation. In MPSK, the binary code word is used 
to vary the phase of the carrier. There are therefore m = 
2**k different phases in MPSK. These sequences of binary 
digits need not be orthogonal as in Orthogonal Binary Phase 
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Shift Keying. The length of the transmitted symbol is only 
k bits long as compared to 2**k bits for Orthogonal BPSK. 
The phases therefore are not all orthogonal leading to 
greater difficulty in decoding and a higher probability that 
an individual code word will be incorrectly decoded, 
especially as m becomes large. 

As in other types of Phase Shift Keying, MPSK 
can be represented as in equation 5.8 where theta = phase = 
(2pi/ni) i^ i*0^ 1, . . .m* 1. 



vc<t) = A cos (wt + theta delta) 



(egn 5.8) 



An advantage of MPSK is that as long as the 
symbol transmission rate does not increase, the carrier 
bandwidth does not increase even as the size of the binary 
code word increases. The disadvantage, as stated earlier, 
has tc dc with error rates in decoding. Decoding requires a 
phase coherent reference of considerable stability 
especially as m increases and respective phases become 
closer together. Practical limits for this type of 
modulation is m = 8 or a 3 bit binary code word due to the 
complexity of the encoding and decoding hardware. 

2. Ampli t ude ^ift Key ing ( ASK ) 

Amplitude Shift Keying, as the name implies is 
modulation of the carrier through variation of its amplitude 
due tc variations in the baseband waveform. The information 
in the baseband waveform is thereby imparted to the carrier 
through this modulation. Since the amplitude of the carrier 
varies between successive symbols or binary code words, the 
power also varies. Recovery is by comparison of the 
transmitted power during each symbol to the possible m = 
2**k power levels. 



a. Multiple Amplitude Shift Keying (MASK) 

MASK is a type of group encoding where k bits 
are combined into a single waveform for transmission and 
subseguent recovery. Since each symbol or binary code word 
contains k bits, there are m = 2^*k different symbols and 
conseguently m different amplitudes of the carrier possible 
during the symbol duration time. An assignment scheme which 
maps the binary code words into the m different amplitudes 
would be to simply space them Ai = A/mi volts apart, where A 
is the maximum amplitude and mi represents the decimal 
equivalent of the binary code word from i = 1 to 2**k. 

The analytic form of a MASK waveform is 
represented in equation 5.9 where all the variables are the 
same a presented earlier except for Ai which is equal to one 
of i different amplitudes depending on the bit sequence to 
be transmitted. 

vc(t) = Ai cos(wt +. delta) (egn 5.9) 

A 

MASK is not very popular in satellite 
communications since it depends on a carrier of very stable 
amplitude. This is not generally practical under actual 
conditions. However, MASK does have the spectral advantage 
of a constant bandwidth of 2/Ts, Ts equal to symbol 
duration, even as the binary code word increases in length. 

b. Quadrature Amplitude Shift Keying (QASK) 

Quadrature Amplitude Shift Keying (QASK) is a 
hybrid digital signal modulaticn technique. It represents 
separate amplitude modulation of the quadrature components 
of a common carrier. In that sense it is both ASK and PSK. 
Implementation of QASK involves simultaneous application of 
M-ary Amplitude Shift Keying to the quadratures. See 
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equation 5.10, where A1 and A2 are derived according to an 
assignment scheme as in MASK. 

vc(t) = A1 cos(wt + delta) + (egn 5.10) 

A2 sin (wt + delta) 



By this method, the effective bit rate is 
doubled over that of ordinary MASK. In other words, QASK 
results in the same data rate as modulating 2k bits onto the 
carrier at the same time or during the same symbol duration 
with MASK. The advantage to QASK lies in the fact that this 
increase in bit rate comes at no further increase in 
bandwidth. The disadvantages associated with MASK are still 
present with an additional increase in the complexity of the 
decoding equipment. Decoding of QASK must be done in two 
steps. First the signal is recovered in phase through a 
phase-coherent correlator and then each amplitude modulated 
signal is recovered as before in MASK by comparison of the 
power levels of the received signal. 

3- Frequency ^ift Keyi ng ( FSK ) 

Frequency Shift Keying (FSK) , is the generic term 
used to describe a number of digital signal modulation 
techniques which cause variations in the carrier frequency 
by interaction with the baseband waveform. Decoding is 
generally through measurement of the frequency of the 
received signal. 

a. Mimimum Shift Keying or Fast Frequency Shift 
Keying 

Power requirement for the transmission or 
retransmission of a satellite signal is extremely important. 
It is desireable to limit power required to accurately 



transmit a piece of information to the minimum possible 



consistent with 


allowable 


error rates. 


For 


this reason. 


detection of FSK 


signals 


is limited to 


coherent methods as 


power required 


increases 


significantly 


for 


non-coherent 


met hods. 

MSK 


or FFSK 


are identical 


and 


represent a 



modulation technique known as continuous phase FSK. They 
get their names from the fact that '’fast” indicates more 
bits per second can be transmitted in a given bandwidth than 
ordinary BPSK and '’minimum" refers to the minimum modulation 
index for which orthogonal signalling occurs [Ref- 7 ]- 

Analytically/ FFSK can be expressed as in 
equation 5.11 and equation 5.12. 

vc(t) = A cos((wc + delta w)t) (eqn 5.11) 

vc(t) = A cos(±deltaw t) cos (wc t) - 5.12) 

A sin (ideltaw t) sin(wc t) 



As is noted/ FFSK can be envisioned as the 
summation cf separately modulated in-phase and quadrature 
components of the carrier by the baseband waveform. The 
technique is similar to that used in OQPSK- The frequency 
variation in the carrier waveform is between wc + deltaw and 
wc - deltaw. Maximum separation in phase of these two 
frequency components occurs at pi as noted in equation 5.13, 
where Tb represents the bit or symbol duration. 

(w2-w1)Tb = pi (eqn 5.13) 

This can be converted to the modulation index 
mentioned above simply by converting to frequency as in 
equation 5.14. 
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h = (f2-f1)Tb = .5 



(eqn 5. 1 4) 



For FFSK, the frequency deviation or separation 
between frequencies of the transmitted carrier is exactly 
l/2Tb and the deviation from the carrier frequency is 1/4Tb. 
FFSK can therefore be rewritten from equation 5.12 as in 
equation 5.15 by substitution for deltaw = 2pi(1/4Tb). 

VC (t) = A cos(±pi t/2Tb) cos (2pi fc t) - (egn 5.15) 

A sin dpi t/2Tb) sin(2pi fc t) 

The error rate in this modulation technique is 
identical to that of BPSK. Frequency Shift Keying can be 
implemented to incorporate any given modulation index by 
simple calculation of the frequency separation with equation 
5.14 and substitution into equation 5.15. For M-ary 
Frequency Shift Keying, the separation between frequencies 
must be at least 1/Ts, where Is is the symbol duration, in 
order to avoid carrier energy from one frequency being 
incorrectly interpreted as carrier energy from another 
frequency. This results in a modulation index of 1. Figure 
5.16 represents the phasor diagram of an FFSK modulated 
signal. 

4. Quadrature pa rtia l Respons e Sig nalli ng (QPRS) 

Quadrature Partial Response Signalling (QPES) 
represents a specific type of the general class of digital 
signal modulation techniques called Partial Response 
Signalling. Partial Response Signalling is a method in 
which a controlled amount of intersymbol interference (ISI) 
is allowed during the encoding process. Previously any 
overlap in successive signals resulted in aliasing and an 
increased probability of error in decoding. The idea was 
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Figure 5. 16 Phaser Diagram of FFSK or MSK 

first introduced by Lender in 1963 [Ref. 12]. His concept 
was that in knowing and controlling the amount of 
interference allowed during the encoding process, 
compensations can then be made for the ISI at the receiver. 
Allowing for a limited amount of ISI makes it possible to 
transmit at rates egual to the Nyquist rate, something that 
is not possible with many other systems [Ref. 13]. 

QPES is implemented, as with previous quadrature 
systems, through simultaneous modulation of the quadrature 
components of a common carrier. However, this time the 
modulation is accomplished with a PRS system. The 
modulation of the respective quadrature components 
represents the impulse response of one of a number of linear 
filters to the bits cf the baseband waveform. The linear 
filter employed depends on the class of the Partial Response 
Signalling system being used- There are 5 classes of linear 
filters commonly used in PES and they are illustrated 
analytically in Table 1 and graphically in Figure 5.17. 

Output of a QPBS system can be represented as in 
equation 5.1b, where an and bn represent the nth bit to be 
modulated and h(t - nis) represents the contribution of the 



79 




1 



TABLE 1 

QPfiS SYSTEM POLYNOMIALS 



SYSTtM POLYNOMIAL 
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Figure 5-17 QPHS Linear Filter Impulse Responses 

impulse response of the linear filter at time t for bit n. 
[Ref. 14] 
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(eqn 5. 1 6) 



oo 

7c(t) = yan h(t - E Ts) cos (wt + delta) 
•'*•00 

^ tn h(t - n Ts) sin(wt + delta) 
n»*oo 



QPRS has the advantage of rapid transmission rates 
with no increase in bandwidth and excellent error handling 
performance. The cost of such performance improvement comes 
in the form of a higher required signal to noise (SNR) ratio 
when compared to other binary systems. 
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VI. COHP OTBB SIMilTION OF D IGI TAL 
SIGJAl HOPaLATlON TECHNIfiOES 

Due to the nature of digital signal modulation, it is 
possible to simulate systems through the use of the digital 
computer. This chapter is a description of the computer 
simulation included in the Appendix of this thesis. 

A. GEBEEAL DESCHIPTICN OF THE PBOGRAM 

The program was constructed using top-down design and 
modular programming, FORTRAN was selected as the 
programming language due to its capanili ties in the area of 
numerical operations and the availability of mathematical 
routines to be used as modules in the program. Testing was 
accomplished on each module as the program was developed. 
The program was designed to be used interactively in order 
to allow for instructional use. Although the program was 
made user friendly in that error checking is accomplished on 
user inputs, care must be taken to insure instructions are 
followed correctly. 

The program consists of twenty modules, sixteen of which . 
were written by the author, three which were taken from the 
Double Precision International Mathematics and Statistics 
Library (IMSLDP) and one which was taken from a NON-IMSL 
library which resides on the IBM 3033 located at the Naval 
Postgraduate School. The IBM 3033 was used exclusively for 
development and testing. 

The development was accomplished using the WATFIV 
compiler; however, the program was written to operate with 
the VS FORTRAN compiler as well. Testing and operation was 
done using the VS FORTRAN compiler. Figure 6.1 is a block 



diagram of the relationship of the modules for a 
representative modulation subroutine. 

B. BlIH CONTBOL I1000I2 

MAIN operates as the control module for the entire 
computer simulation and accomplishes limited output. It 
introduces the program and allows the user to input various 
parameters of the digital signal modulation technique to be 
simulated and various other control functions. These inputs 
inc 1 u de : 

1. Modulation technique 

2. Class of QPRS system (when appropriate) 

3. Digital logic scheme 

4. Baud or symbol rate 

5. Bits/binary code word (when appropriate) 

6. Carrier frequency 

7. Number of samples to be generated 

8. Carrier max amplitude 

9. Initial phase angle 

10. Ose of random number generator or no 

11. Seed for RNG (when appropriate) 

12. Number of repetitions of the simulation 

Once the user has input the desired characteristics of 
the modulation to be simulated, MAIN calls the appropriate 
module to begin the actual calculation. The values 
necessary to calculate the time series of the signal are 
passed to the subroutines as parameters. MAIN calls the 
required subroutine the number of times the user specifies 
as repetitions. Each call to a subroutine produces the 
Discrete Fourier Transform and amplitude spectrum of the 
signal. It is these respective amplitude spectrums which 
produce the statistics in the final output after MAIN calls 
the statistics subroutine the final time. 
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Figure 6,1 Representative Block Diagram 
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C. 50BBOOTINE BPSK 

This description of SUBEOOTINE BPSK, a module that 
simulates Binary Phase Shift Keying, will serve as a guide 
on the construction and operation of subsequent modules 
involved in the actual signal modulation. A complete 
description including its interaction with other program 
subroutines will be included- Since other signal modulation 
subroutines interact in the same manner within the program, 
subsequent descriptions will concentrate on differences to 
this basic module. 

SDBROUTINE BPSK begins with variable declarations as 
with all FORTEAN programs. All variables are declared to be 
either type integer or double precision. All calculations 
in the program are carried out in double precision accuracy- 

variable initialization is next- Significant variable 
initialized at this point are TIME, the time of the first 
sample which is set to 0- STEP, the delta time between 

samples is set to the normalized Nyquist sampling rate- 

OMEGA, the angular frequency and DELTA, the initial phase 
angle are computed in radians- NBAUD, a variable which 

keeps track of the number of symbols which have been 

modulated, is set to 1 or the first symbol is being 

modulated- NBAUD remains at its .present value until total 
elapsed time exceed 1 symbol duration or BAUDD, the 

normalized baud rate. 

Subsequently , an initial call is made to SUBROUTINE 

STREAM to receive the value of the first bit to be 

modulated- The modulation is then carried out according to 
equation 6 - 1 - 

The values of each sample are assigned to an array for 
storage along with the corresponding time increment- 

Modulation continues with the first drawn bit until 1 bit or 
symbol duration is exceeded when STREAM is called to draw 
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vc<t) = A v(t) cos ( Ht + delta) 



(egn 6.1) 



A = amplitude 

\r(t) = bit to be modulated 
w = angular frequency 
t = time 

delta = initial phase angle 

another bit. This entire process is repeated until the 
number of samples specified by the user in MAIN is reached. 

At this point, only on the first repetition of the 
simulation, SUBROUTINE PLOT is called which gives the user 
the opportunity to plot the time series of the simulated 
signal. Upon return from PLOT the subroutine calls the IMSL 
routine FFTEC to generate the Discrete Fourier Transform of 
the time series produced. This is accomplished on each 
repetition of the program. 

Again on the first repetition only, information on the 
number of the principle harmonic of the FFT is displayed. 
Plot is called again to give the user the opportunity to 
plot the amplitude spectrum. The values of the FFT and the 
amplitude spectrum are computed on each successive 

repetition of the program and the amplitude spectrum is 
added to the statistics being accumulated by a call to 
SUBROUTINE STATS. 

D. SUBROUTINE DBPSK 

SUBROUTINE DBPSK simulated Differential Binary Phase 
Shift Keying and is essentially the same as SUBROUTINE BPSK. 
The difference lies in the fact that what is encoded during 
the modulation process is the product of the random bit 
drawn and the value of the bit previously modulated. A 
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reference bit equal to +1 is used to determine the first bit 
to be modulated. In this manner, modulation of •►I means no 
change occurs between successive bits and the phase of the 
carrier remains the same. Modulation of a -1 indicates a 
change of bits has occurred and is indicated by a change of 
phase of the carrier between successive bit durations. The 
rest of the module remains unchanged. 

E. SOBBOOTIHE OBPSK 

Orthogonal Binary Phase Shift Keying is accomplished in 
this module. Again SUBROUTINE OBPSK is constructed 

essentially the same as SUBROUTINE BPSK. In Orthogonal BPSK 
the only bit streams that are modulated are n bits long 
representing each binary code word, where n = (k = 

number of bits per word is limited to 6 bits in MAIN). Each 
sequence of n bits is orthogonal to every other allowed 
sequence. 

The way this is accomplished in this module is first to 
draw a series of k random bits from SUBROUTINE STREAM. This 
series of k 1*s and 0*s is then changed to its decimal 
equivalent from 0 to 2**k-1 and saved for later use. 

The program continues with a call to SUBROUTINE ORTHO 
which generates an n X n orthogonal matrix of +1*s and -1*s 
through the use of the Hadamard matrix and the Kronecker 
product [fief. 15]. The Hadamard matrix is a 2 X 2 
orthogonal matrix of 1's and - 1 ' s. The Kronecker product is 
the matrix which is formed when a matrix is expanded to 
twice its original size by replacing each element of the 
Hadamard matrix by the product of the Hadamard matrix and 
the original matrix. The value of the decimal equivalent to 
the k random bits is used to identify the row of the matrix 
and the respective columns represent the value of the binary 
digit which is modulated onto the carrier by BPSK. This 
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sequence continues until either one bit duration is 
exceeded, when another orthogonal bit is received from ORTHO 
or when all the orthogonal bits in a row are modulated- Then 
another sequence of k random bits is converted to its 
decimal equivalent and a new row and column of the 
orthogonal matrix is identified. Execution continues in 
this manner until all required samples have been produced- 
It should be noted that the bit rate for Orthogonal BPSK is 
n times the baud rate. This significantly increases the 
signal bandwidth requirement and decreases the time 
increment between successive samples in the simulation- The 
rest of the module is the same as those previously 
present ed- 

F. SOBfiCOTINE QPSK 

Quadrature Phase Shift Keying, simulated in SUBROUTINE 
QPSK, is accomplished through simultaneous BPSK modulation 
of the quadrature components of the carrier. Two random 
bits are drawn from successive calls of SUBROUTINE STREAM 
and each bit is in turn modulated onto the carrier 
components and the sum formed. This process is repeated at 
time intervals equal to the normalized Nyquist sampling rate 
until TIME exceeds one symbol duration, BAUDD. At this time 
two new random bits are drawn and the modulation continues 
until all required samples are produced. The rest of the 
modules is the same as those before. Equation 6-2 
represents the analytical expression used to simulate 

G- SOBROOTINE OQPSK 

SUBROUTINE OQPSK is essentially the same as SUBROUTINE 
QPSK with these minor deviations- The time cf the first 
sample is artificially initialized as TIME = -5(BAUDD) or 
1/2 the symbol duration- In this manner, the two random 
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vc(t) = A v1(t) cos(wt + delta) + 
A v2(t) sin (wt delta) 



(egn 6.2) 



v1 (t) = in-phase bit 
v2(t) = quadrature tit 



bits first being modulated are both known to have existed on 
the carrier for 1/2 the symbol duration. The first of thest 
two bits only remains on the carrier until TIME = 1 (BAUDD) 
when a third random bit is drawn and modulation begins with 
this bit. The second random bit is allowed to remain until 
TIME = 1.5 (BAUDD) when a fourth random bit is drawn to 
replace it and so on. In this manner the two bits being 
modulated are offset in the time they change by 1/2 the 
symbol duration. This necessitates keeping track of the 
number of bits modulated on both the in-phase and quadrature 
component of the carrier. The remaining portion of the 
module is the same as SUBROUTINE QPSK. 

H. SUBROUTINE MPSK 

This module simulates M-ary Phase Shift Keying. This is 
accomplished by changing the phase of the carrier to any one 
of n =2**k phases where k is the number of bits in the 
binary code word- The max length of the binary code word is 
limited to 10 in the MAIN program. This modulation is 
accomplished by use of equation 6.3 to simulate the MPSK 
system. 

The program first draws k random bits and converts them 
to the decimal equivalent. At this point modulation begins 
at TIME = 0 according to equation 6.3. Modulation continues 
until one symbol duration is exceeded when a new set of k 
random bits is drawn, conversion to decimal takes place and 
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vc(t) = A cos[ wt + (2pi m) /n ♦ delta] 



(egn 6.3) 



m = decimal equivalent of binary code word 
n = 2**k; k = bits in binary code word 



modulation continues or until all desired samples have been 
produced. The program continues as in previous subroutines. 

I. SOBBODTIBE HASK 

The principles cf operation of SUBROUTINE MASK is 
similar to SUBROUTINE MPSK which also involves a group 
encoding system. SUBROUTINE MASK allows for M-ary Amplitude 
Shift Keying. The analytic expression used in the 

simulation is provided in equation 6.4 Ai is one of 2**k 
equally spaced amplitudes. 

VC (t) = Ai cos(wt + delta) (eqii 6.4) 

Once again a set of k random bits is drawn . from 
SUBROUTINE STREAM. The decimal equivalent of the k bits is 
computed and used to adjust the amplitude Ai according to 
the assignment routine expressed in equation 6.5. 

Ai = A/m (egn 6.5) 

A = max amplitude 

m=1 ton; n=2**k 



This process is repeated at each symbol duration until 
all desired samples have been produced. The rest of the 
module remains unchanged. 
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J. SOBBOOTINE QASK 



SOBROUTINE QASK combines the elements of the previously 
described amplitude modulation technique with the now 

familiar quadrature modulation technique. Two separate 

streams of random numbers are generated and converted to 
their decimal form. These numbers are used to produce two 
separate amplitudes with the same assignment scheme used in 
MASK. These amplitudes modulate the quadrature components 
of a common- carrier and the output is formed by the sum of 
the quadratures. Modulation continues in the above manner 
until a symbol duration elapses when two new amplitudes are 
calculated. Processing terminates when all required samples 
are computed. 

K- SUBROUTINE MSK 

Minimum Shift Keying is frequency shift keying in which 
the modulation index is 1/2. The modulation index is 
represented in equation 6.6 

h = .5 = (deltaf) Tb (eqn 6.6) 

deltaf = frequency separation 

Tb = tit duration 



The frequency deviation from the carrier is .5 (deltaf) = 
1/4Tb. Converting this frequency deviation to radians 
yields pi/2Tb. This angular frequency deviation from the 
central carrier is either ±pi/1Tb depending on the value of 
the bit to be modulated. 

The simulation in SUBROUTINE MSK is accomplished by 
drawing a random bit and generating the signal according to 
equation 6-7. 



VC (t) = A cos[ (w ±pi/2Tb) t + delta] 



(e^n 6.7) 



Once again modulation continues until one bit duration 
is elapsed when another bit is drawn. Modulation stops when 
all required samples are produced. The rest of the 
subroutine remains unchanged. 

L. SOBBOOTINE HFSK 

This subroutine modulates a signal simulating M-ary 
Frequency Shift Keying. The separation between adjacent 
frequencies is established as 1/Ts, where Ts is the symbol 
duration. This amounts to a modulation index equal to 1 . 
Initially the entire range of frequency deviation from the 
carrier is calculated as (n-l)Ts. The mean frequency 
deviation is then calculated. In other words the range of 
frequencies of the modulated signal is the frequency of the 
central carrier ±1/2 the magnitude of the entire range of 
frequency deviation. The minimum frequency is used as the 
base for computing the frequency to be used to modulate the 
signal during each respective symbol duration. 

At this point in the subroutine a series of k random 
bits <k<10) is drawn and converted to decimal. This decimal 
number is multiplied by the frequency separation and added 
to the minimum frequency and converted to radians. The 
analytical expression of the signal produced with this 
subroutine is illustrated in equation 6.8. 

vc(t) = A cos[ (w + 2pi mf) t * delta] (egn 6.8) 

mf = modulation frequency 



This modulation continues until a symbol duration 
elapses when a new frequency deviation or mf is calculated 
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or mcdalation stops do to completion of all required 
samples. 

H. SOBBOaTINE QPBS 

This module simulates 5 different classes of Quadrature 
Partial Response Signalling systems. The technique employed 
involves the summation, at each time increment, of all 
responses to the linear filter representing the class. This 
is accomplished for all of the bits being modulated during 
the duration of the simulation. 

Initally, two arrays of random bits are generated 
representing those bits which would be modulated onto the 
quadrature components of the carrier during the length of 
time the signal is to be simulated. The the impulse 
response to each of these n bits is calculated for the time 
of the respective sample. These impulse responses represent 
the modulation for the bit in question onto the respective 
portions of the carrier and the sum is formed according to 
equation 6.9. 

oo 

VC (t) = h (t - nTs) cos (wt + delta) + (eqn 6.9) 

CO ns-oo 

y A h (t - nTs) sin (wt + delta) 
n=-co 

h (t - nTs) = impulse response of nth bit at time t 



The time of the first sample is initialized at 6(BAUDD) 
in order to ensure that when the first sample is taken, 
contributions from bits modulated at time 0 are also summed 
do to the overlapping nature of the QPR signals. The value 
of the last sample is also formed by the sum of the 
quadrature components existing at 6(BAUDD) after it is 
produced. The rest of the program operation remains the 
same as previous modules. 
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H. SOBfiOOTIHE STREAH 






This subroutine interacts with the IMSL random number 
generator SUBROUTINE GGUBFS to produce a random number 
between 0 and 1 and make assignment on the basis of the 
value of this random number to random bits according to the 
digital logic scheme specified in the MAIN program. It also 
enables the user to manually insert binary digits if that 
option was specified previously in MAIN. 

O. SUBROUTINE ORTHO 

CBTHO produces an n X n orthogonal matrix from the 
Hadamard matrix by forming successive Kronecker products 
[Ref. 15]. It also selects the appropriate row and column 
of the orthogonal bit to be modulated and passes this bit 
back to SUBROUTINE OEPSX. Point of entry to SUBROUTINE 

ORTHO is controlled by a flag set in OBPSK. 

P. SUBROUTINE PLOT 

PLOI is an interactive module that allows the user to 
determine whether or not a graph of the output of program 
calculations is to be produced. The user is also able to 
selectively vary certain parameters associated with the 
plot. PLOT calls SUBROUTINE DTPLOT which actually produces 
the graph and performs the output. UTPLOT is a NON-IMSL 
library routine from the Naval Postgraduate School. 

Q. SUBROUTINE STATS 

The statistical calculations associated with the 
amplitude spectrums of the various functions generated on 
successive repetitions of the program are computed in 
SUBROUTINE STATS. Upon completion of each repetition, each 
module calls STATS where the sum, sum of the squares, sum of 
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the cubes and sua of the guartics of each element of the 

amplitude spectrum of the signal produced are computed and 

stored. When the last repetition of the program is 

complete, dilH instructs STkTS to compute and output the 

final statistics, i.e., mean, variance, skewness and 
kurtosis of each component of the amplitude spectrum 
according to the estimations contained in [fief. 16 ]. 
Additionally the mean variance and variance of the variance 
is calculated and output. Point of entrance to this 
subroutine is controlled by flags set in the individual 
Modulation subroutines or in : 1 AIN. 
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VII. S TATISTI CS OF THE FAST FOORIER TRA N SFO RMS 



A. DESCRIPTION OF THE PROBLEM 

As Stated in Chapter II, if it can be shown that the 
statistics of the FF T can be somehow linked to a particular 
digital signal modulation technique, then the modulation 
technique can be identified on the basis of those statistics 
alone. Although it is beyond the scope of this thesis to 
actually derive those relationships if they exist, it seem 
prudent to attempt to determine if there is statistical 

differences between the components of the FFT*s as they are 

derived in the enclosed computer simulations. 

The statistic chosen to do the hypothesis testing is the 
F-test since the true mean and variance of the distribution 
need not be known [Ref. 17]. In addition, the information 
necessary to calculate the F-statistic is readily 
accumulated in SOBROOTINE STATS. Calculation of the 

F-statistic necessitated the writing of a small computer 
program to be used for that purpose once the output from the 
main program was generated and reformatted. 

B. AHALYSIS-OF-VARI ANCE 

If there were no differences attributable to the 

modulation technique, then it would be reasonable to assume 
that for a certain set of characteristics, the components of 
successive FFT*s would be the same. In other vords, if a 
signal was modulated by BPSK and a statistically significant 
number of FFT*s were generated of the signal, then the 
statistics of those FFT*s would be assumed to be related. 
On the other hand, if it can be shown that the statistics 
are not related (i.e., not from the same distribution), then 
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significance can be placed in the variation among modulation 

m 

techniques. 

1 , The F-Di s tributi on 

The assumptions necessary to use the F-distribution 

are : 

1 . Normal distribution 

2. Random samples 

3. Independent samples 

These assumptions are not too difficult to intuitively 
accept in the model that will be proposed. Once again it 
would be expected that the FFI*s of successive identically 
modulated signals to be related. For the simulations 
conducted as part of this test, they were all generated 
using bits from a random number generator for which it can 
be shewn that each sequence of bits passes statistical tests 
for randomness and independence. The assumption of 
normality could also be argued on the basis of the central 
limit theorem. 

The value of the F-statistic with which comparison 
will be made is 1.51 for 14 degrees of freedom in the 
numerator, an infinite number of degrees of freedom in the 
denominator and a 90^ confidence region- How these values 
were obtained will be detailed under the design of the 
experiment. 

2- Desig n of ;^e Sxperime nt 

The experimental data used in the computation of the 
statistics was chosen to minimize the random contributions 
of variables other than the modulation technique. The 
computer simulations included in this thesis were used to 
generate the statistics and a separate program also included 
in the appendix was used to calculate the F-statistic. The 
variables in the experiment include: 
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1 . Modulation technique 
« 

2. Logic type 

3. Baud rate 

4. Bits per binary code word 

5. Carrier frequency 

6. Carrier amplitude 

7. Initial phase angle 

8. Number of samples generated 

9. Time between samples 

10. Seed for random number generator 

11. Number of repetitions or trials 

Fifteen different modulation techniques were 
compared. In all cases bipolar logic was simulated, the 
baud rate was held constant at 1200 baud and the maximum 
carrier amplitude was established at 1 volt. Additionally 
the phase angle of all simulations was 0 degrees and the 
seed for the random number generator was 1 thereby ensuring 
the same random sequence of bits. The number of bits per 
binary code word did vary among the modulation techniques. 
When the modulation technique was M-ary, the bits per code 
word was always 3 so it is possible to infer that bits per 
code word is a function of the modulation technique. The 

time between samples was the normalized Nyguist sampling 
rate. This did vary between modulation techniques but was 
necessary to derive an accurate FFT. Carrier frequency also 
varied between simulations but was always twice the lowest 
carrier frequency recommended in the computer simulation. 
In most cases this was 2400 Hz, or twice the baud rate. 
Finally each simulation was repeated 100 times and the 
statistics of the FFT*s based on those 100 repetitions. The 
number of samples produced was always 6.4 so there are 33 
components of the FFT. 

What this amounts to can be illustrated through an 
Analysis-of-Variance table as shown in Table 2 [Eef. 17]. 
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TABLE 2 

AHALYSIS-OF-YARIANCE TABLE 



1 



n Observations in Each of r Groups 







Observations 


Sum 


Mean 




1 






Yix 


Yi„ 


^i. 


Yi. 




■*11 


* 1 A 


■•u ..... 


* i/i 




2 


I'll 


Yiz 


I'll 


. Y^n 


Yz. 


Yz. 


Group 


3 


>'31 




>'33 


■ Y^ 


Yy 


Yy 




1 

1 

1 

1 










r 


Yrl 


Yr2 




. Yrr, 


Yr 


Yr 






* n 




rj ..... 


m 


' r. 


■* r. 



In th€ case of the simulation described above, this amounts 
to 100 observations (repetitions or trials) from each of 15 
groups (modulation techniques) for each of the 33 components 
of the FFT. The degrees of freedom in the numerator is 
equal to (15-1) or 14 while the degrees of freedom in the 
denominator is equal to 15(100-1) or 1485. Table values for 
the F-d istribution use infinity as the degree of freedom for 
values greater than 120. 

An F-test was also performed to determine if there 
is a relationship among the mean variance, mean skewness or 
mean kurtosis of the 15 different modulation techniques. 
Again the assumption is made that these statistics would be 
normally distributed from modulation techniques which had 
the same or statistically similar components of their FFT*s- 
In this case there were again 15 modulation techniques to 
compare (14 degrees of freedom) but only 33 observations 
(the variance, skewness or kurtosis of the respective 
elements of the FFT). This number of observations yields 
15(33-1) = 480 degrees of freedom in the denominator. 
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3 . Hypothesi s a d d Hypothesis Test i ng 

The hypothesis posed for the model is that the means 
of the individual components of the FFT’s from the 15 
modulation techniques are from the same distribution. 
Therefore^ the respective means must be equal. Also tested 
is that the mean variance, mean skewness and mean kurtosis 
of all components of the FFT’s of each modulation technique 
are equal. These hypothesis are illustrated in equations 
7.1 through equation 7.4. Each hypothesis is tested and 
compared to the F-distribution selecting a rejection region 
or level of significance of .9. This equates to an 
F-statistic of 1.51 for 14 degrees of freedom in the 
numerator and an infinite number of degrees of freedom in 
the denominator. 

4 . Results 

The means of the 33 components of the FFT for each 
of the 15 modulation techniques were compared using th^e 
F-distribution. The results are shown in Table 3, Table 4 
shows the F-statistics associated with the comparison of the 
mean of the variance, mean of the skewness and mean of the 
kurtosis for the respective groups. In addition a complete 
summation of the results of the F-test are included in the 
appendix. 

C. CCHCIOSIOH 

The results of the F-test indicate that the null 
hypothesis must be rejected at the .9 significance level and 
the alternate hypothesis accepted that the means are not 
equal for FFT components 9-14. In addition the F-statistic 
for the comparison of mean variance, mean skewness and mean 
kurtosis all fall in the rejection area. 
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Ho; XBAR11 = XBAfi2 1 = • • • = XBARij (egn 7.1) 

Hi: XBAR11 * XBAR2 1 • • • # XBARij 

i = ith modulation technique 
j = jth component of the FFT 

Ho: MVAR1 = MVAR2 = • • • = MVAR15 (eqn 7.2) 

HI: MVAR1 # MVAR2 ^ aVAR 15 

HO: MSKEW1 = aSKEW2 = • • • = MSKEW15 (eqn 7.3) 

HI: aSKEWI HSKEW2 ••• 4 MSKEW15 

Ho: MKOR1 = MK0E2 = • • • = HKORIS (eqn 7.4) 

HI: MKUE1 # MKUE2 ^ • 4 BKOEIS 



This indicates that there is a difference between the 
statistics of of the FFT*s of the respective modulation 
technique. Since all the parameters used in the generation 

of these statistics were held essentially constant among the 
trials^ it can be assumed that these differences are due to 
the modulation technique employed. It has not been 

established yet what those differences may be. This is an 
area where future research and study is warranted. 

The results of this test would seem to point to those 
components of the FFT which offer the best opportunity to 
develop those relationships or differences. FFT components 
9-14 have obvious differences; however, FFT components 20, 
24 and 25 also have large F-statistics but do not fall in 
the rejection area. These components may also be 
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TABLE 3 

P-STATISTICS OF FFT COMPONENTS 



FFT 

CCMPCNENT 


AMONG 
G£OUP SUM 


WITHIN 
GEOUP SUM 


TOTAL 
OF SQUARES 


F-STAT 


1 


0.121 


17.69 


17.81 


0.727 


2 


0.112 


18.48 


18.59 


0.642 


3 


0. 108 


20.65 


20.76 


0.557 


4 


0. 1 35 


25.40 


25.54 


0.562 


5 


0. 172 


32.74 


32.91 


0.558 


6 


0.254 


41.72 


41.98 


0.647 


7 


0.422 


50-10 


50.52 


0.893 


8 


0.962 


82.83 


83.79 


1.231 


9 


7.591 


390.2 


397.8 


2.063 


10 


49.09 


1681.0 


1730. 1 


3.097 


11 


53. 01 


1732. 1 


1785. 1 


3.246 


12 


61.56 


2071.7 


2133.3 


3.152 


13 


45. 73 


1832.8 


1878.6 


2.646 


14 


3.864 


250.5 


254.4 


1.636 


15 


0.539 


98.13 


98.67 


0-582 


16 


0.366 


110-1 


110-5 


0.353 


17 


0.5 04 


15 1-0 


151. 5 


0.354 


• 18 


0.943 


140.3 


141.2 


0.713 


19 


1.299 


158.5 


159.8 


0.370 


20 


2.67 1 


20 0. 7 


203. 4 


1.412 


21 


1.853 


226.3 


228.2 


0.868 


22 


2.277 


232. 1 


234.3 


1.041 


23 


1.7 12 


197.6 


199.3 


0.919 


24 


2.5 12 


18 1.9 


184.4 


1.476 


. 25 


2. 1 44 


17 1.3 


173.5 


1-327 


26 


1.257 


143.7 


145-0 


0. 928 


27 


1.233 


113.9 


115. 1 


1.149 
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TABLE 3 

F-STATISTICS OF FFT COilPONENTS (con*t) 



28 


0.858 


11 5.3 


116.7 


0.785 


29 


0.547 


63.72 


64.27 


0.911 


30 


0.356 


55.77 


56.13 


0.677 


31 


0.183 


41.41 


51.59 


0.469 


32 


0.051 


13.79 


13.84 


0.389 


33 


0.031 


10.92 


10.96 


0.3 34 



TABLE 4 

F-STATISTICS OF THE MEAN 7A2IANCES, 
MEAN SKEWNESS AND MEAN K3RT0SIS 





AMONG 


WITHIN 


TOTAL 






GROOP 


SUM 


GROUP SUM 


OF SQUARES 


F-SIAT 


VARIANCE 


3. 144 


E5 


3.755 E6 


4.070 E6 


2. 370 


SKEWNESS 


5. 307 


E10 


7.031 Ell 


7.562 Ell 


2.588 


KURTOSIS 


8. 280 


Ell 


1. 262 E13 


1.344 E13 


2.250 



interesting to examine. In addition, it should be noted 
that the statistics of the FFT associated with each 
respective modulation technique also display some striking 
differences. These statistics may prove in some way to be a 
fingerprint of the modulation techniques themselves. 
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appendix a 

FORTRAN PROGRAM FOR DIGITAL SIGNAL MODULATION 
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MODULATION TECHNIQUE = ORTHOGONAL BPSK 
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APPBHDII E 
BESOLTS OF F-TEST 
F-STATISTICS FOR IHE 1 COMPONENT 
SOM OF SOUS = 0.2091487882294994D 02 
SUM OF SUMS**2 = 0 . 18 1 06554 69453404D 02 
SOM OF SQUARE OF SOMS = 0. 4 12809751 0679048D 02 
AMONG GROOP SUM OF SQUARES = 0. 1 21 1883 1361545 15D CO 
WITHIN GROUP SOM OF SQUARES = 0, 1769 3744 943 466 1 4D 02 
TOTAL SUM OF SQUARES = 0 . 17 81 4933257 08 159D 02 
DEGREES OF FREEDOM IN NUMERATOR = 0, 1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0- 148500 000 COOOOOOD 04 
F-STATISTIC = 0.72650632 14387600D 00 

F-STATISTICS^ FOR THE 2 COMPONENT 

\ 

SUM OF SUMS = 0.2305792717805020D 02 

SOM OF SUMS**2 = 0- 18947.0123401529 ID 02 

SUM OF SQUARE OF SUMS = 0. 46638251 26 655468D 02 

AMONG GROUP SUM OF SQUARES = 0- 1 1 1 93 71 755000 360D 00 

WITHIN GROUP SOM OF SQUARES = 0. 184806298274 8736D 02 

TOTAL SUM OF SQUARES = 0 . 1 85925670 02 98739D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 400000000000000D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 1485000000000000D 04 

F-STATISIIC = 0.64247464650 14737D 00 



F-STATISTICS FOR THE 3 COMPONENT 

SOM OF SUMS = 0.2693734535827635D 02 

SUM OF SUMS**2 = 0 . 21 246 74 9 21 974 1 86D 02 

SUM OF SQUARE OF SUMS = 0. 5921566381 358867D 02 

AMCNG GROUP SUM OF SQUARES = 0, 1 08 40 9 588 1685 1 83D 00 

WITHIN GROUP SUM OF SQUARES = 0 . 20 65 4 592581 6 0598D 02 

TOTAL SUM OF SQUARES = 0 . 20 76300216977449D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 400000 0 00000000 D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 1 4850000000000 OOD 04 

F-STATISTIC = 0.5567362242775778D 00 
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F-STATISTICS FOR THE 4 COHPONENT 

SUH OF SUMS = 0-3 175688631627065D 02 

SOM OF SOMS**2 = 0. 262 1 1 045 4230 021 4D 02 

SUM OF SQUARE OF SUMS = 0. 8 0705 1 3093 726266D 02 

AMONG GROUP SUM OF SQUARES = 0. 1 34 71 80 903696 0 15D 00 

WITHIN GROUP SOM OF SQUARES = 0. 25403994 11362951D 02 

TOTAL SUM OF SQUARES = 0. 255387122039991 ID 02 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 400000000000000D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 1 4850000000000 OOD 04 

F-STATISTIC = 0-5624997485041869D 00 

F-STATISTICS FOR THE 5 COMPONENT 

SUM OF SUMS = 0-3669104349531 116D 02 

SUM OF SUMS**2 = 0- 33806337 403 17695D 02 

SUM OF SQUARE OF SUMS = 0- 1 069682473356158D 03 

AMCNG GROUP SUM OF SQUARES = 0. 1721940248396 148D 00 

WITHIN GROUP SUM OE SQUARES = 0- 32 73 66 549298 2079D 02 

TOTAL SUM OF SQUARES = 0 . 32 9088489 5466 04 ID 02 

DEGREES OF FREEDOM IN NUMERATOR = 0- 1 4000000 000 OOOOOD 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0- 1485000000000000D 04 

F-STATISTIC = 0-5579331866788867D 00 

F-STATISTICS FOR THE 6 COMPONENT 

SUM OF SUMS = 0-4 1 60097878.94843 ID 02 

SUM OF SUMS**2 = 0 - 43 1 32 366 42 02 0094D 02 

SUM CF SQUARE OF SUMS = 0- 1 40 82 1 53 58 007360D 03 

AMCNG GROUP SUM OF SQUARES = 0- 2544544005 1 1 9 440D 00 

WITHIN GROUP SUM OF SQUARES = 0. 41 724 1 51 062 1 9358D 02 

TCTAL SUM OF SQUARES = 0 - 4 1 97860 546270552D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0, 1 4000000 OOOOOOOOD 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0- 1485000000000000D 04 

F-STATISTIC = 0-646875756 1623442D 00 
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F-STATISTICS FOfi THE 7 COBPONENT 

SOW CF SUMS = 0.4810387092348 11 ID 02 

SOM OF SOMS**2 = 0 . 52064 093 8 16 1 1 74 4D 02 

SOM OF SQOARE OF SOMS = 0. 1 964653659903027D 03 

AMONG GROOP SOM OF SQOARES = 0. 4219987280210724D 00 

WITHIN GROOP SOM OF SQOARES = 0.50099440 1562 1441D 02 

TOTAL SOM OF SQOARES = 0 . 50 52 143888 42354 8D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 4 00 0000 000 00000 D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 1485000000000000D 04 

F-STATISTIC = 0.8934632362547178D 00 

F-STATISTICS FOR THE 8 COMPONENT 

SUM OF SUMS = 0. 5.623255941290002D 02 

SUM CF SUMS**2 = 0. 8590 2856 22098 20 3D 02 

SUM OF SQUARE OF SUMS = 0.3070018365343552D 03 

AMONG GROUP SUM OF SQUARES = 0. 96 1 95 1 2065933 320D 00 

WITHIN GROUP SUM OF SQUARES = 0. 8283283785563848D 02 

TOTAL SOM OF SQOARES = 0 . 83 7947 89062 23 1 8 1 D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 400000 0000 00 OOOD 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 1 4 8500 000 00000 OOD 04 

F-STATISTIC = 0. 1 2318247369 14027D 01 

F-STATISTICS FOR THE 9 COMPONENT 

SUM OF SUMS - 0.9 931605295569187D 02 

SUM OF SUMS**2 = 0 .40 4353 10 56 1 7588 1 D 03 

SUM OF SQUARE OF SUMS = 0, 1 4 1 664 24 1 728 31 50D 04 

AMCNG GROUP SUM OF SQUARES = 0. 7590638589699638D 01 

WITHIN GROUP SUM OF SQUARES = 0 . 39 0 1 8668 1 444 7566 D 03 

TOTAL SUM OF SQUARES = 0 . 39 77773200344562D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 400000000000000D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0, 1 48500 000 00000 OOD 04 

F-STATISTIC = 0.2063499133280512D 01 
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F-STAIISTICS FOB THE 10 COMPONENT 

SOM OF SUMS = 0.2029942628781924D 03 

SOM OF SOMS**2 = 0. 1757564757819924D 04 

SOM OF SQOAfiE OF SOMS = 0. 7 65569225045 1 960D 04 

AMONG GROOP SOM OF SQOARES = 0. 4908580866354581D 02 

WITHIN GEOOP SOM OF SQOARES = 0. 16 8 1 00 78353 1 540 4D 04 

TOTAl SOM OF SQOARES = 0 . 173009 3643978950D 04 

DEGREES OF FREEDOM IN NOMERATOE = 0 . 1 4000000 000 00000 D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 1485000000000000D 04 

F-STATISTIC = 0. 30S73096841 92641D 01 

F-STATISTICS FOE THE 11 COMPONENT 

SOM OF SOMS = 0. 2004773302247246D 03 

SOM CF SOMS**2 = 0. 18 1 18 750 7756 1 6 1 7D 04 

SOM OF SQOAfiE OF SUMS = 0- 7 9804398631 08266D 04 

AMCNG GEOOP SOM OF SQOARES = 0- 5301 029200839379D 02 

WITHIN GROOP SOM OF SQOARES = 0. 1732070678930534D 04 

TOTAL SOM OF SQOARES = 0 . 1 785080970938928D 04 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 4000000 OOOOOOOOD 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 1485000000000000D 04 

F-STATISTIC = 0.3246332537532912D 01 

F-SIATISTICS FOR THE 12 COMPONENT 

SUM OF SUMS = 0, 22 10834652773109D 03 

SUM CF SUMS**2 = 0 . 21 658597 652 1 894 4D 04 

SUM OF SQUARE OF SUMS = 0, 9 414279932 820438D 04 

AMCNG GROUP SOM OF SQUARES = 0. 61557 53358218842D 02 

WITHIN GROOP SUM OF SQUARES = 0. 2071 71 6965890740D 04 

TOTAL SUM OF SQUARES = 0 - 21 3327449 94 72928D 04 

DEGREES Of FREEDOM IN NUMERATOR = 0. 1 400000000000000D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 148500000 OOOOOOOD 04 

F-STATISTIC = 0. 3 151731454585565D 01 
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F-STATISTICS FOR THE 13 COMPONENT 

SUM OF SUMS = 0- 1 979583434471791D 03 

SOM OF SUMS**2 = 0. 190470 1075347160D 04 

SOM Of SQUARE OF SUMS = 0. 7 185446672 840862D 04 

AMCNG GROUP SOM OF SQUARES = 0. 4 57 29 46 290 1 50 773D 02 

WITHIN GROUP SUM OF SQUARES = 0. 18 328466086 1 875 ID 04 

TOTAL SUM OF SQUARES = 0 . 18 7857607 1 520259D 04 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 400000000000000D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 1485000000000000D 04 

F-STATISTIC = 0.2646478671459854D 01 

F-STATISTICS FOR THE 14 COMPONENT 

SUM OF SUMS = 0- 9 2277362726 14619D 02 

SUM OF SUMS*=«'2 = 0.26008271 S6886254D‘ 03 

SUM OF SQUARE OF SUMS = 0. 9 5407268 1 908 1 072D 03 

AMCNG GROUP SUM OF SQUARES = 0. 386398570461 9236D 01 

WITHIN GROUP SUM OF SQUARES = 0. 2505 4 1 992869 544 4D 03 

TOTAL SUM OF SQUARES = 0 . 254405978574 1 636D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 400000000000000D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0 . 1 48500 000 COOOO OOD 04 

F-STATISTIC = 0. 1 6358873774 98235D 01 

F-SIATISTICS FOR THE 15 COMPONENT 

SOM OF SUMS = 0,68809263362 10145D 02 

SUM OF SUMS**2 = 0. 10 18300362534556D 03 

SUM Of SQUARE OF SUMS = 0.3 69506 28634 15098D 03 

AMCNG GROUP SUM OF SQUARES = 0. 538 58 63 804 584 070D 00 

WITHIN GROUP SUM OF SQUARES = 0. 98 13497339004056D 02 

TOTAL SUM OF SQUARES = 0 . 9867355977049896D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0- 1 40000000000000GD 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 1485000000000000D 04 

F-STATISTIC = 0.5821433971075596D 00 
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F-5TATISTICS F0£ THE 16 COMPONENT 

SOM OF SUMS = 0,69613873264504310 02 

SUM OF SUMS**2 = 0. 1 1 371 5703396946 3D 03 

SUM OF SQUARE OF SUMS = 0. 3 597053597 878584D 03 

AMCNG GROUP SUM OF SQUARES = 0. 3663260306209 388D 00 

WITHIN GROUP SUM OF SQUARES = 0. 1 10 1 1 8649799 0677D 03 

TOTAL SUM OF SQUARES = 0. 1 1 04849758296886D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 400000000000000D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 1485000000000000D 04 

F-STATISTIC = 0, 3528623485827816D 00 

F-STATISTICS FOR THE 17 COMPONENT 

SOM OF SUMS = 0.7052854165691613D 02 

SUM OF SUMS**2 = 0. 1 5480420049 14454D 03 

SUM OF SQUARE OF SUMS = 0;. 38203241 48469269D 03 

AMCNG GROUP SUM OF SQUARES = 0. 504 1 4 06896350 342D 00 

WITHIN GROUP SUM OF SQUARES = 0. 150983876342976 ID 03 

TOTAL SUM OF SQUARES = 0. 15 14880170326 1 1 1D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 400 0000 000 00 000 D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 1 48500 000 00000 OOD 04 

F-SIATISTIC = 0.3541763825767676D 00 

F-STATISTICS FOE THE 18 COMPONENT 

SUM OF SUMS = 0. 79504528933 26030D 02 

SUM OF SUMS*=f'2 = 0 . 1454 1 276 56 12229 ID 03 

SUM OF SQUARE OF SUMS = 0- 5 157 1 61626920556D 03 

AMCNG GROUP SUM OF SQUARES = 0. 9 43 1 8 1 546320 8 067 D 00 

WITHIN GROUP SOM OF SQUARES = 0- 1402556039853086D 03 

TOTAL SUM OF SQUARES = 0 . 14 1 1 98 78553 16294D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 4000000 OOOOOOOOD 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0- 1 485 00 000 00000 00 D 04 

F-STATISTIC = 0.7 1330207S7581566D 00 
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F-STATISTICS FOR THE 19 COHPONENT 

SUM OF SUMS = 0.8206200282896035D 02 

SUM OF SUMS**2 = 0. 1642503769158388D 03 

SUM OF SQUARE OF SUMS = 0-5788746875066006D 03 

AMONG GROUP SUM OF SQUARES = 0. 1 299298669532 475D 01 

WIIHIN GROUP SUM OF SQUARES = 0. 15846 1 6300407728D 03 

TOTAL SUM OF SQUARES = 0. 15976092871 03052D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 4000000 000 OOOOOD 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 14850000000000 OOD 04 

F-STATISTIC = 0.8697276809711277D 00 

F-STATISTICS FOR THE 20 COMPONENT 

SUM OF SUMS = 0.9520048302882186D 02 

SUM OF SUMS**2 = 0 . 20940 865 584 1 4699D 03 

SUM OF SQUARE OF SUMS = 0 . 8 71321 92 1 88 1 5827D 03 

AMONG GROUP SUM OF SQUARES = 0. 267 1 1 31 239535 162D 01 

WITHIN GROUP SUM OF SQUARES = 0. 20 069 5436622 654 1 D 03 

TOTAL SUM OF SQUARES = 0 . 20 3366567 86 2 1 893D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 400000000000000D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 1 485000000000000D 04 

F-STATISTIC = 0, 1 411744637781583D 01 

F-STATISTICS FOR THE 21 COMPONENT 

SUM OF SUMS = 0.89941399521 50146D 02 

SUM of' SUMS**2 = 0.2335862791765991D 03 

SUM OF SQUARE OF SUMS = 0 . 7 2457556355 1 0632D 03 

AMONG GROUP SUM OF SQUARES = 0- 1 85278 5403586 404D 01 

WITHIN GROUP SUM OF SQUARES = 0. 2263 40 52354 1 0884D 03 

TOTAL SOM OF SQUARES = 0 . 22 81 933089446748D 03 

DEGREES Of FREEDOM IN NUMERATOR = 0. 1 4 000000 000 00 000 D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0 . 1 48500 000 C 0000 OOD 04 

F-STATISTIC = 0. 8682828488687504D 00 
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F-STATISTICS FOE THE 22 COMPONENT 

SUM OF SUMS = 0. 1 035648300735959D 03 

SOM OF SUMS**2 = 0- 24148 1 0868436007D 03 

SUM OF SQUAEE OF SUMS = 0. 9 42734 1493765320D 03 

AMONG GROUP SUM OF SQUARES = 0. 227689214 1650 1 19D 01 

WITHIN GROUP SOM OF SQUARES = 0. 2320 537453498354D 03 

TOTAL SUM OF SQUARES = 0 . 23 4330 637 49 14855D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1400000000000000D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 148500000 OOOOOOOD 04 

F-STATISTIC = 0. 1 040764077320930D 01 

F-STATISTICS FOR THE 23 COMPONENT 

SUM OF SUMS = 0.9031234717154746D 02 

SUM OF SUMS**2 = 0. 2047768 1 16906234D 03 

SUM OF SQUARE OF SUMS = 0. 7 14937975597 1323D 03 

AMCNG GROUP SOM OF SQUARES = 0- 1 711833054881 912D 01 

KITHIN GROUP SUM OF SQUARES = 0. 1976 27 43 1 934 652 1 D 03 

TOTAL SUM OF SQUARES = 0- 1993392649895340D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 400 00 00 000 00 OOOD 02. 

DEGREES OF FREEDOM IN DENOMINATOR = 0- 1485000000000000D 04 

F-STATISTIC = 0.9 1878 22552243536D 00 

F-STATISTICS FOR THE 24 COMPONENT 

SUM OF SUMS = 0- 9 171706807335522D 02 

SUM OF SUMS**2 = 0. 1900438745495707D 03 

SUM CF SQUARE OF SUMS = 0- 8 1395624 803 4 03 89D 03 

AMCNG GROUP SUM OF SQUARES = 0. 253 1 548763025 405D 01 

WITHIN GROUP SUM OF SQUARES = 0 . 1 8 1 9 04 31 206 9 230 3D 03 

TOTAL SUM OF SQUARES = 0 . 1 8 4435860 832255 7D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 4 00 0000 000 00 000 D 02 

DEGREES CF FREEDOM IN DENOMINATOR = 0- 1 48500000 OOOOOOOD 04 

F-STATISTIC = 0. 1476188171339999D 01 
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F-STITISTICS FOfi THE 25 COMPONENT 

SOM OF SUMS = 0.8440126623444366D 02 

SUM OF SUMS**2 = 0- 17822 548 95437147D 03 

SUM OF SQUARE OF SUMS = 0. 689290 1073945840D 03 

AMONG GROUP SUM OF SQUARES = 0. 2 1 4385 1 9 1 2627 546D 01 

WITHIN GROUP SOM OF SQUARES = 0. 171332 58 846 97689D 03 

TOTAL SUM OF SQUARES = 0 . 1734764403823964D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 400 0000 00000 000 D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 1485000000000000D 04 

F-STATISTIC = 0- 13272514998 40134D 01 

F-STATISTICS FOE THE 26 COMPONENT 

SUM OF SUMS = 0-6923579472492147D 02 

■SUM OF SUMS**2 = 0 . 148 1 9 850 5836 6 125D 03 

SUM OF SQUARE OF SUMS = 0 . 4 452980334 1859 1 3D 03 

AMCNG GROUP SUM OF SQUARES = 0. 1 257250 1 5339 1 605D 01 

WITHIN GROUP SUM OF SQUARES = 0. 1437455255024266D 03 

TOTAL SUM OF SQUARES = 0 . 1450027756558 182D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 400 0000 000 00 000 D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 1485000000000000D , 04 

F-STATISTIC = 0.92773892874 7691 OD 00 

F-STATISTICS FOfi THE 27 COMPONENT 

SUM OF SUMS = 0.6366092822129744D 02 

SUM OF SUMS**2 = 0 . 1 1 7790 68 1 99 1 273 OD 03 

SUM OF SQUARE OF SUMS = 0 .3 S351590 41535181D 03 

AMCNG GROUP SUM OF SQUARES = 0. 1 2333 49 853537 058D 01 

'killHIN GROUP SUM OF SQUARES = 0. 11 3855 522949 7378D 03 

TOTAL SUM OF SQUARES = 0 . 1 1 50888728032 74 9D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 400000000000000D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0 . 14850 0000 00000 OOD 04 

F-STATISTIC = 0. 1 149027974258136D 01 
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F-STATISTICS FOR THE 28 COMPONENT 

SUM OF SUMS = 0.5387592494868525D 02 

SOM OF SOMS**2 = 0 . 11 8 6 1 334 42790 129D 03 

SOM OF SQUARE OF SUMS = 0. 27925959579 1 5495D 03 

AMCNG GROUP SUM OF SQUARES = 0. 85751 9098531 2 5 17D 00 

WITHIN GROUP SUM OF SQUARES = 0. 1 1582074 832 1 0974D 03 

TOTAL SOM OF SQUARES = 0 . 1 1 667826741 96286D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 400000000000000D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 14850000000000 OOD 04 

F-STATISTIC = 0.7853366268738306D 00 

F-STATISTICS FOR THE 29 COMPONENT 

SUM OF SUMS = 0.4536144151906700D 02 

SOM OF SUMS**2 = 0.65643051 39058856D 02 

SUM OF SQUARE OF SUMS = 0. 1 91 90 53830677748D 03 

AMONG GROUP SUM OF SQUARES = 0. 5472802462192574D 00 

WITHIN GROUP SUM OF SQUARES = 0. 63723 9975599 1 08 1 D 02 

TOTAL SUM OF SQUARES = 0 . 64 27 12778 06 13 007D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 400 0000 000 00000 D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0- 1 48500 000 00000 OOD 04 

F-S.TATISTIC = 0.9 1097231 44851782D 00 

F-STATISTICS FOE THE 30 COMPONENT 

SOM OF SUMS = 0.3 732568029970994D 02 

SUM OF SUMS**2 = 0. 570576522863441 9D 02 

SUM OF SQUARE OF SUMS = 0 . 1 284820507967475D 03 

AMCNG GROUP SUM OF SQUARES = 0. 3560 162347433721D 00 

WITHIN GROUP SUM OF SQUARES = 0 . 5577 283 1 7783 767 2D 02 

TOTAL SUM OF SQUARES = 0 . 56 1288480 1 3 1 2 009D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 400000 OOOOOOOOOD 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 1485000000000000D 04 

F-STATISTIC = 0,6 770886363437515D 00 
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F-STAIISTICS FOR THE 31 COMPONENT 

SUM OF SUMS = 0-2973876286269883D 02 

SUM OF SUMS**2 = 0 - 42 178 336 1065906 8D 02 

SUM OF SQUARE OF SUMS = 0. 77262052049 18280D 02 

AMONG GROUP SUM OF SQUARES = 0- 18302 45094226 048D 00 

WITHIN GROUP SUM OF SQUARES = 0. 41 40 57155860 9885D 02 

TOIAl SUM OF SQUARES = 0 - 4 1 58874 0095 52 1 45D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 400 0000 000 00 OOOD 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0, 1485000000000000D 04 

F-STATISTIC = 0- 4688645251806333D 00 

F-STATISTICS FOR THE 32 COMPONENT 

SUM OF SUMS = 0. 21 42367689889270D 02 

SUM OF SUMS**2 = 0 . 1 4 1 48 204 254 1 6 03 8D 02 

SUM OF SQUARE OF SUMS = 0- 3 56 57234 02 1 828 90D 02 

AMONG GROUP SUM OF SQUARES = 0. 505 897 1 897285 643D- 0 1 

WITHIN GROUP SUM OF SQUARES = 0. 1379 1631 9 139 420SD 02 

TOTAL SUM OF SQUARES = 0 - 1 3 84222 1 6329 1 4 95D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0, 1 400000000000000D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0- 1485000000000000D 04 

F-STAIISTIC = 0.3890854828465455D 00 

F-STATISTICS FOR THE 33 COMPONENT 

SUM OF SUMS = 0. 1 7'32004830070321D 02 

SUM OF SUMS*^2 = 0 . 1 1 1 55 820 9706648 ID 02 

SUM OF SQUARE OF SUMS = 0- 231 2899538600082D 02 

AMCNG GROUP SUM OF SQUARES = 0. 3 130057 176754676D-01 

WITHIN GROUP SUM OF SQUARES = 0. 10 92453 1 0 168 0480D 02 

TOTAL SUM OF SQUARES = 0 - 1 0 95583 1 58857235D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 4000000 OOOOOOOOD 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 1485000000000000D 04 

F-STATISTIC = 0.3039120267386334D 00 
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F-STATISTICS OF MEAN VARIANCE, SKEWNESS AND KOETOSIS 

F-STATISTICS FOR THE 1 COMPONENT 

SUM OF SUMS = 0- 1 1 116174394 22250D 05 

SUM OF SUMS**2 = 0- 431 9425777749392D 07 

SUM OF SQUARE OF SUMS = 0. 1 861260733 960752D 08 

AMONG GROUP SUM OF SQUARES = 0, 3 143833877400 904D 06 

WITHIN GROUP SUM OF SQUARES = 0. 375540737351 886 1 D 07 

TOTAL SUM OF SQUARES = 0 . 40 6979076 1258 952D 07 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 4 000000 000 00000 D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 480000 OOOOOOOOOOD 03 

F-STATISTIC = 0.28702236365 19026D 01 

F-STATISTICS FOR THE 2 COMPONENT 

SUM OF SUMS =-0. 3062348822555926D 07 

SUM OF SUMS**2 = 0 .7751 235739069674D 12 

SUM OF SQUARE OF SUMS = 0. 2376364 1 81 020636D 13 

AMONG GROUP SUM OF SQUARES = 0. 5306 562 1 0 1 878 76 1 D 11 

WITHIN GROUP SOM OF SQUARES = 0. 703 1 12538 1 1 8 4634D 12 

TOTAL SOM OF SQUARES = 0 - 756 178 1 59 1 372509D 12 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1 4000000000 OOOOOD 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0, 4800000000000000D 03 

F-STATISTIC = 0. 25876266202 17669D 01 

F-STATISTICS FOR THE 3 COMPONENT 

SUM OF SUMS = 0. 1 087692759375071D 08 

SUM OF SUMS**2 = 0 . 1 36 82 3 55 71 001 09 OD 14 

SUM OF SQUARE OF SUMS = 0. 35212589426623 16D 14 

AMONG GROUP SUM OF SQUARES = 0, 8 28 04 30 04 8704 6 94D 12 

WITHIN GROUP SUM OF SQUARES = 0. 1 261 530754556 778 D 14 

TOTAL SUM OF SQUARES = 0 , 13 44335055043825D 14 

DEGREES OF FREEDOM IN NUMERATOR = 0- 1 400 0000 000 00 000 D 02 

DEGREES OF FREEDOM IN DENOMINATOR = 0. 48000000000000 OOD 03 

F-STATISTIC = 0.22504442145S6077D 01 
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